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A computer-aided data acquisition system was developed 
and a microthermocouple probe constructed to obtain thermal 
distributions in turbulent buoyant jets exposed to a cross- 
flowing ambient fluid. The system performed high speed 
temperature measurements as a microthermocouple probe was 
automatically traversed through a sequence of preprogrammed 
positions under the control of a microcomputer. Operability 
of the apparatus was demonstrated by measuring temperature 
distributions in planes perpendicular to the streamwise axis 
of jets from which contour plots of temperature were gener- 
ated. Using temperature distributions along with velocity 
distributions allow buoyant jet characteristics to be 
computed, including the entrainment rate of ambient fluid, 
jet trajectory, and heat transfer to the ambient. The 
experimental technique is discussed and temperature contour 
plots for a jet at various planes are presented. 
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INTRODUCTION 



I . 

Buoyant jets are very common in nature. We see them in 
the form of exhaust gases emitted from smoke stacks of 
refineries, mills and ships. We see them in the form of 
heated waste water expelled into the sea from power plants 
and from the main propulsion condensers in steam driven 
ships and submarines. It is no wonder that the fluid 
mechanics and heat transfer characteristics of buoyant jets 
have been of interest to environmental, civil and mechanical 
engineers for decades. To evaluate their ecological impact, 
and of most recent interest, to harness buoyant jets as a 
means of detecting military targets and guiding weapons, it 
is necessary to develop models which accurately predict 
their trajectory and decay. 

Most studies to date have dealt with buoyant jets rising 
through a quiescent ambient fluid; however, in nature most 
problems involve flowing ambient fluids. Relatively little 
experimental work has been done with buoyant jets in cross- 
flow and, according to Hilder [Ref. 1], the trajectories of 
jets and the entrainment rates of ambient fluid predicted 
from previous work do not agree well with one another. Most 
mathematical models of buoyant jets in crossflow assume 
Gaussian profiles for velocity and temperature. Nickodem 
[Ref. 2] has shown through experiments that in fact, the 
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Gaussian profiles of velocity are altered by crossflow. 

This leads one to suspect that the same may be true for the 
temperature profiles. 

The objective of this work was to develop a system to 
thermally map a buoyant jet in crossflow. Then, by measur- 
ing both velocity and temperature distributions, improved 
computations of entrainment, trajectory aad heat transfer 
characteristics of jets can be made thereby giving rise to 
more accurate models . 
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II. BUOYANT JETS DISCHARGED TO A CROSSFLOW 



A. PROPERTIES OF BUOYANT JETS 

A buoyant jet is characterized by a momentum and a 
density differential between the jet and its surrounding 
ambient resulting from a variation in temperature and/or 
fluid concentrations. Therefore, fluid motion in the jet 
is governed by both inertial and buoyant forces. The non- 
dimensional ratio of these forces, known as the densiometric 
Froude number, provides an important quantitative measure- 
ment of jet characteristics and is shown below. 

U2 

^ " gD(Pa -Po^/Po 

where U^ is the jet's discharge velocity, g is the accelera- 
tion of gravity, D is the discharge diameter of the jet, 
is the density of the crossflowing ambient and is 
the density of the jet fluid at its point of discharge. 

The Gaussian velocity and temperature profiles assumed 
by most models of buoyant jets are very similar. Velocity 
behavior is given by: 

U(r) = Uj^ exp(-r^/B^) 

where U is the centerline velocity, r is the independent 
m 

variable and a radial distance from the centerline of the 
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jet, and B is defined as nominal jet halfwidth. As r 
approaches B, velocity decays to (l/e)U^ [Ref. 3]. Simi- 
larly, temperature behavior is given by: 

T(r) = exp (-r^/A^B^) 

where is the centerline temperature, r and B are defined 
the same as above and A, a spreading ratio, is the inverse 
of the turbulent Schmidt number (S) . S is defined as the 
ratio of the molecular momentum and mass dif f usivities and 
is equal to where v is the kinematic viscosity and 

D is the binary mass diffusion coefficient associated with 
substances A and B [Ref. 4]. Although A varies inversely 
with the Froude number, the change is very slight, and in 
the case where substances A and B are both water, A is 
slightly greater than 1. Hirst [Ref. 3] found A to vary 
between 1.16 at F = 0 to 1.11 at F = infinity. The net 
effect then, is a more gradual temperature decay than was 
found with velocity. 

Most buoyant jet models consider the entrainment of the 
ambient fluid into the jet and are based on relevant conser- 
vation equations of mass, momentum and energy. In conserva- 
tion of mass, the downstream change in total mass of the jet 
is equated to the mass of the entrained fluid. The conser- 
vation of momentum must consider both vertical and horizontal 
contributions. Changes in vertical momentum are equated 
to the buoyant forces while changes in the horizontal 
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momentum of the jet are equated to the horizontal momentum 
of the entrained fluid. The conservation of energy involves 
energy changes resulting from variations in the ambient 
temperature as caused by the jet. Hilder iRef. 1] developed 
the following governing equations in non-dimensional differ- 
ential form. 



CONTINUITY 



HORIZONTAL MOMENTUM 



VERTICAL MOMENTUM 



ENERGY 



(u b ) 
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B . FLOW REGIMES 

The jet passes through several regimes as it travels 
from the nozzle through the ambient. The three regions most 
frequently referred to are shown in Figure 1. They are the 
zone of flow establishment, the zone of established flow 
and the far-field zone [Ref. 3]. In the zone of flow estab- 
lishment, the velocity and turbulei^ce profiles transform from 
the conditions within the nozzle to a free turbulent flow 
condition. It is in this region that the jet begins to mix 
with the ambient fluid; however, the flow is still more 
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Figure 1. Typical Flow Regions in a Buoyant Jet 



strongly influenced by the nozzle discharge conditions than 
by the ambient. When the turbulent mixing has reached the 
centerline of the jet, the zone of established flow is said 
to begin. In this region, the profiles have assumed their 
free turbulent shapes. Now the flow is governed by the 
jets' momentum and buoyancy as well as by the condition of 
the crossflow. The far field zone is defined as that region 
in which jet momentum is depleted and the jet fluid is con- 
vected and diffused by the ambient currents and turbulence. 

C. EFFECTS OF CROSSFLOW 

At the immediate exit of a cylindrical nozzle, a verti- 
cally discharged buoyant jet has a nearly uniform velocity 
distribution and has the same cross-sectional shape as the 
nozzle itself. The velocity gradient between the jet and 
the crossflowing ambient creates longitudinal shear stresses 
at the jet's sides, a positive pressure region immediately 
upstream and a negative pressure region immediately downstream 
of the jet. This results in the deflection of the jet's 
trajectory in the downstream direction (Figure 2), the 
creation of counterrotating vortices at the jet's outer 
edges and the deformation of the original circular cross- 
sectional shape into the form of a kidney. As the stream- 
wise axis of the jet approaches the direction of the cross- 
flow, these effects become progressively less pronounced. 
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Figure 2. Coordinate System for a Buoyant Jet 




III. EXPERIMENTAL APPARATUS 



A. SYSTEM OVERVIEW 

A surplus milling machine was configured with synchron- 
ous drive motors interfaced with a microcomputer that auto- 
matically positioned its bed. It was used as a three-dimensional 
positioning platform in the same manner as in the laser 
Doppler velocimetry work undertaken by Nickodem [Ref. 2] . 

The milling machine was placed adjacent to a rectangular 
plexiglass flume through which the crossflowing ambient 
fluid flowed. A vertical nozzle was installed in the base 
of the flume to provide the jet. A temperature probe was 
suspended through an opening in the top of the flume above 
the nozzle by an arm attached to a base mounted on the 
milling machine bed as shown in Figure 3. As the probe was 
automatically traversed through a series of preprogrammed 
positions across the jet, temperature data was automatically 
sensed and stored at high speeds by the computer. 

B. CROSSFLOW SYSTEM 

As illustrated in Figure 4, the crossflow circulation 
pump took water from the cylindrical 248.8 1 (65.7 gal) 
reservoir shown in Figure 5 and discharged through 5.076 cm 
(2 in) diameter tubing into a cylindrical flow settling 
chamber 30.46 cm (12 in) in diameter and 60.91 cm (24 in) 
tall located within the 60.91 cm x 60. 91 cm x 88.83 cm 
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Figure 3. Probe and Traversing Mechanism 
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Figure 4. Crossflow and Jet Loop Piping Diagram 
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Figure 5. Flume Circulation and Figure 6. Inlet and Flow 

Jet Loop Pumps Settling Chambers 




(24 in X 24 in x 35 in) inlet chamber shown in Figure 6 . The 
settling chamber was sealed at its bottom so that the water 
spilled from its top into the inlet chamber through honey- 
combed flow straighteners to reduce turbulence and evenly 
disperse the flow. To further reduce turbulence, the flow 
was broken by another stack of honeycombed flow straighteners 
and a layer of fiberglass filter material located immediately 
above the normal operating water level. The flow next 
entered a 24.4 cm x 32. 39 cm x 182.9 cm (9.625 in x 12. 75 in x 72 in) 
flume shown in Figure 7 through a vertical section of the 
same honeycombed material mentioned above. To avoid inadver- 
tent spillage over the sides of the flume during system 
start-up, a 5.076 cm (2 in) diameter overflow pipe was 
located in the inlet chamber. During normal operation, 
a gate valve in this piping was closed. The flow left the 
flume through a 7.614 cm (3 in) diameter pipe at its end 
and re-entered the crossflow circulation pump reservoir. A 
gate valve located in this piping and shown in Figure 8 was 
used to regulate the water level and flow velocity in the 
flume. The optimum adjustment of this valve was determined 
by trial and error to be closed two turns from its fully 
open position. Either a globe or ball valve would have been 
more appropriate for this purpose; however, neither was 
readily available, so the gate valve was used. The bracket 
shown at the base of the flume in Figure 8 maintained align- 
ment between the flume and the milling machine. The water 
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Figure 7. Flume Arrangement 



Figure 8. Flume Discharge Piping 
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in the flume was gradually heated by repetitive circulation 
through the crossflow pump and by the addition of the heated 
water from the jet. To maintain a constant temperature 
crossflow, cooling water from a refrigerated bath shown in 
Figure 9 was circulated through a coil of 1.269 cm (.5 in) 
diameter copper tubing located in the flow settling chamber. 
Also, fresh water was added at the flow settling chamber 
as an equal amount was drained from the crossflow pump 
reservoir through a 1.269 cm (.5 in) diameter pipe. Cross- 
flow temperature was monitored by a Type-T thermocouple 
located in the inlet chamber. 

C. JET SYSTEM 

In reference to Figures 4 and 5, the jet flow pump 
circulated water from a rectangular 26.27 1 (6.94 gal) 
reservoir and discharged through 1.26 cm (.5 in) diameter 
tubing to a 33.0 cm x50.8 cmx54.6 cm (13 in x20 inx21.5 in) 
head tank (Figure 10) . The amount of flow to the head tank 
was regulated by a globe valve. Due to a low flow rate 
to the head tank, water was also recirculated back to the 
reservoir in order to maintain sufficient flow through the 
jet pump to prevent overheating it. A constant water level 
was maintained in the head tank by a stand pipe which allowed 
overflow back to the reservoir. Sufficient flow into the 
tank was maintained to make sure that it slightly overflowed 
continuously. Water drained from the bottom of the head 
tank through 1.26 cm (.5 in) diameter tubing and passed 
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Figure 9. Refrigerated Bath 
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Figure 10 . Head Tank 
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through a rotometer, a .95 cm (.375 in) tubing reducer, a 
dye injection system (Figure 11) , a water heater (Figure 12) 
consisting of approximately 6.09 m (20 ft) of .95 cm (.375 
in) diameter copper tubing coiled in a heated bath and 
finally a 7.144 cm (.28125 in) nozzle which discharged into 
the bottom of the flume. Flow was controlled by pinching 
the tubing between the heater and the nozzle with surgical 
clamps. Drainage from the crossflow reservoir discussed in 
Section III.B was used to replenish the jet reservoir. The 
dye injection system, used in photographing the jet, was 
located approximately 8.23 m (27 ft) upstream of the nozzle 
to minimize any disturbance to the jet that it might have 
caused. The majority of this distance was taken up by the 
heating coil mentioned above. The vertical distance between 
the top of the stand pipe in the head tank and the tip of 
the nozzle in the flume was 2.2 m (86.5 in) which equated to 
21.56 KPa (3.127 psig) . Jet flow temperature was monitored 
by a Type-T thermocouple located within the jet flow tubing 
approximately 1.167 m (46 in) from the nozzle. 

D. TEMPERATURE PROBE 

Measuring temperatures in a buoyant jet with a thermo- 
couple is intrusive. To reduce the probability of distorting 
results, steps were taken to minimize the cross-sectional 
area of the temperature measuring device as seen by the flow 
of the jet. A .0254 mm (.001 in) diameter Type-E micro- 
thermocouple was selected. The suspension device for the 
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Figure 11 . Dye Injection System Figure 12 . Jet Heater 








microthermocouple had to be rigid and have a small cross- 
sectional area, for reasons discussed above, as well as be 
an electrical insulator to prevent interference with the 
thermocouple performance. A glass annulus approximately 
1.45 mm (.057 in) in diameter and 11.27 cm (4.4375 in) in 
length was chosen. One lead of the thermocouple was 
threaded through the annulus and the other was glued with a 
fast drying modelers' glue along the outer surface, allowing 
the microthermocouple junction to protrude slightly from the 
tip of the annulus. The leads at the opposite end of the 
annulus were welded to .0762 mm (.003 in) diameter wire 
which subsequently was connected to 28 AWG extension wire 
to the computer. The annulus was mounted as shown in Figure 
13. Henceforth, this device will be referred to as the 
probe . 

E. PROBE ACTUATOR. ASSEMBLY 

The cross sectional area of the probe as seen by the jet 
was further reduced by orienting the probe tangentially to 
the trajectory of the jet as shown in Figure 14 . This 
photograph indicated that the probe created no noticeable 
interference with the jet hydrodynamics. Probe orientation 
was accomplished by the linkage assembly shown in Figure 15. 
The fixed end of the probe was hinged to a streamlined tube 
23.495 cm (.9.25 in) long with a maximum width and depth, as 
seen by the jet, of 3.175 mm (.125 in) and 6.35 mm (.25 in) 
respectively. It was rigidly connected to the mounting 
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Figure 13. Probe Profile 




Figure 14. Probe in Jet 
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Figure 15. Probe Assembly Linkage 
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bracket as shown in Figure 15 which was connected to the 
probe arm shown in Figure 3 by a single stud which allowed 
pivoting of the probe from side- to-side . The stud was also 
hollowed so that a linkage rod could extend from the hinge 
assembly through the tube and stud to a stroking pulley 
which was rotated by a small motor. The hinge assembly and 
the stroking pulley were spring loaded to reduce hysteresis. 

As shown in Figures 16 and 17, the 1.5 VDC motor, geared to 
one rpm, was directly coupled to a potentiometer as well as 
the drive pulley. The potentiometer was configured in a 
voltage divider such that the amount of motor rotation, and 
ultimately the degree of probe deflection, was proportional 
to the potential difference sensed across the potentiometer. 
Limit switches were installed at the stroke pulley as shown 
in Figures 16 and 18 to prevent damage to the linkage 
assembly due to over-rotation. 

F. MICROCOMPUTER INTERFACE 

The data collection process consisted of adjusting the 
probe angle of deflection, traversing the three-dimensional 
positioning platform and measuring temperature profiles. 

All of the mechanisms which controlled these events were 
interfaced to an HP-9826 computer shown in Figure 19 through 
an HP-6942A multiprogrammer which performed high speed analog- 
to-digital conversions and ultimately provided control 
signals to govern relays within the system. Refer to MAIN_T 
in Appendix B for the microcomputer software which directed 
this process. 
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Figure 16. Probe Actuator Assembly 





Figure 17. Probe Actuator Motor- 

Potentiometer Arrangement 




Figure 18. Probe Actuator 
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Figure 19 . 



HP-9826 Microcomputer 
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1 . Probe Angle Adjustment 



As discussed in Section III.E, a potentiometer con- 
figured as a voltage divider provided probe angle feedback 
to the computer as shown in Figure 20. The direction of 
motor rotation was controlled by the computer through Type- 
MRRlCDL replays connected as shown in Figure 21. When the 
probe was at a desired position, a 5.0 VDC signal was applied 
to pins 6L and 6 r which allowed both relays to assume the 
normally closed (NC) positions which opened the power circuit 
to the motor. When it was desired to rotate the motor clock- 
wise, pin 6R was grounded which resulted in 5.0 VDC applied 
across the coil in the right-hand relay. This caused the 
relay to assiime its normally open (NO) position resulting in 
a 1.5 VDC signal at terminal B of the motor causing it to 
rotate in the clockwise direction. The left-hand relay was 
activated in a similar manner for counterclockwise rotation. 
Refer to PROBE_SUBS in Appendix B for probe positioning 
software . 

2 . 3-D Positioning Platform Movement 

Positioning platform movement was controlled in a 
manner similar to the probe and was discussed in detail by 
Nickodem [Ref. 2] . Refer to MTR_SUBS in Appendix B for 
associated HP-9826 software. 

3 . Temperature Data Collection 

Three thermocouples were monitored in the data col- 
lection process. A Type-T thermocouple located in the inlet 
chamber measured the ambient fluid temperature in the flume. 
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Figure 20. Computer Bus and Probe Voltage Divider Wiring Diagram 
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Figure 21 . Probe Motor Wiring Diagram 



a Type-T thermcouple located in the tubing between the heater 
and the nozzle measured nozzle temperature and a Type-E 
microthermocouple in the probe measured the temperature in 
the jet. The EMF's generated by these thermocouples were 
amplified by "Omega Omni-Amp IIB" millivolt amplifiers shown 
in Figure 22 prior to entering the multiprogrammer for 
analog- to-digital conversion and eventual transformation to 
temperature readings. Fourth-order least squares coeffi- 
cients for this conversion were taken from Beckwith [Ref. 5]. 
Operation of the crossflow circulation pump created suffi- 
cient electrical interference to distort the thermocouple 
signals. This problem was corrected by applying a thin 
coating of silicon sealant to the Type-T thermocouple junc- 
tions and by connecting the crossflow circulation pump 
casing, the nozzle and the jet tubing in the vicinity of the 
nozzle thermocouple to a common ground. Because the jet 
tubing was plastic, it was necessary to manufacture a brass 
"T" connector as shown in Figure 23 which was grounded and 
located in close proximity to the thermocouple junction. 

Refer to T_SUBS in Appendix B for associated software. 
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Figure 22. 



Thermocouple Amplifiers and HP-69 42A 
Multiprogrammer 
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Figure 23. Nozzle Thermocouple Grounding Arrangement 
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IV. EXPERIMENTAL PROCEDURES 



A. CALIBRATION 

Before the data collection process could begin, the 
rotometer, the thermocouples, the probe and the positioning 
platform had to be calibrated. The detailed steps taken 
are discussed below. 

1 . Rotometer 

With a constant level maintained in the head tank 
and the jet tubing disconnected from the nozzle and elevated 
to the same height as the top of the nozzle, five 100 ml 
samples were drawn through the rotometer and timed to the 
nearest 0.01 second at each rotometer reading from 10% 
to 75% in 5% increments. Flowrates and standard deviations 
in ml/s are shown in Table 1. 

2 . Thermocouples 

Since nozzle and probe temperatures were to be 
normalized by the ambient temperature, the nozzle and probe 
thermocouples were calibrated relative to the ambient thermo- 
couple by using the microcomputer program T_CAL in Appendix 

B. The procedure followed is outlined in the initial com- 
ments of the program. Coefficients for first order curve 
fits were solved by the least squares method with the 
mainframe programs TCAL and TFIT found in Appendix C. 
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3. Probe 



As the probe assembly was being developed, it was 
convenient to test its suitability with the probe calibration 
panel shown in Figure 24. Resistance changes across the 
potentiometer were recorded for varying degrees of deflec- 
tion. Analysis of this information led to improved designs 
from the standpoint of reduced hysteresis and repeatability. 
The microcomputer program PROBE_CAL in Appendix B was 
developed to enable calibration of the final design after 
it was installed in the system as shown in Figures 25 and 
26. The calibration procedure is outlined in the preliminary 
comments of the program. 

4 . 3-D Positioning Platform 

The positioning platform was calibrated in a manner 
that placed the tip of the probe at desired locations within 
the flume relative to the tip of the nozzle. Referring to 
the coordinate system illustrated in Figure 2, the center 
of the nozzle was defined as (0,0,0) in xy z-coordinates . 

The following relationships apply to the probe geometry 
shown in Figure 27 : 

X(real) = X - R cos a tan(ir/4 - a/2) 

O 3. 

Y(real) = Y + r (1 - cos y) + R cos a 
op a 

Z(real) = Z - r sin y 
o p 
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Figure 26. Probe Calibration Panel 
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Figure 27. Probe Positioning Geometry 
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(X ,Y ,Z ) was the position of the milling machine bed 
when the tip of the probe was at position (0,0,0) with 
Y = 0 and a = B = 90 degrees. For calibration, the probe 
and probe arm were configured with these settings as shown 
in Figure 28. With measured values of r and R entered 
into the microcomputer program MAIN_T, the calibration was 
accomplished by the program MOTOR_CAL. The step-by-step 
procedure followed was outlined in the subprogram "SUB 
Calibrate." As illustrated in Figures 3, 27 and 28, the 
length of R could be modified to compensate for adjustments 
of a and B to positions other than 90 degrees. Decreasing 
a increased the distance along the Y-axis in which the 
probe could be positioned. Increments of a and B were 
scribed on the top of the base and at the tip of the probe 
arm in Figure 28 to accommodate this change, if desired. 

The program MAIN_T queried the user for the value of a and 
assumed a = B* The calibration software also established 
position limits to prevent driving the probe into the sides 
of the flume. 

B . PRELIMINARIES 

Crossflow velocity was determined by injecting blue food 
coloring into the flow and timing its travel through a 1.0 m 
interval. The average of several trials indicated the 
velocity was .130 m/s (.427 ft/sec) with the flume outlet 
valve closed two turns from its fully open position. 
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Figure 28. Probe Arm Positioning for Motor Calibration 




Figure 29 . 



Typical Buoyant Jet as Obsered with Dye 
In j ected 
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Photographs of jet profiles as shown in Figure 29 were 
taken to determine jet trajectory and halfwidths along the 
streamwise axis. This was done by injecting blue food 
coloring into the jet flow as discussed in Section III.C. 

The specific gravity of the food coloring was found to be 
considerably less than that of water. To eliminate the 
added buoyant effect this would have had on the jet, a small 
quantity of alcohol was mixed with the food coloring as 
suggested by Merzkirch [Ref. 6]. The amount of alcohol added 
was determined by trial and error. As small quantities were 
added and mixed, samples were gently placed on the surface 
of a beaker of water. Pure food coloring laid on the sur- 
face and very slowly mixed with the water. As alcohol was 
added, this buoyant effect grew progressively less and the 
mixture would settle into the water. The mixture was con- 
sidered satisfactory when it no longer laid on the surface, 
but settled to some equilibrium position in the beaker. 

Slide photographs of the jets were projected onto large 
sheets of 3.175 mm (.125 in) grid graph paper and digitized 
along approximate streamwise axes and half-width trajectories. 
A scaling factor was determined by equating the projected 
width of the nozzle to its known outer diameter of 7.9375 mm 
(.3125 in). The above data was fit to the Michael is-Men ter 
Equation [Ref. 7] shown below by the least squares method 
with the mainframe program JETCURV in Appendix C: 
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Correlation coefficients close to 1 were consistently 
obtained. To determine positions within the jet at which 
to make temperature measurements, five evenly spaced posi- 
tions per jet flow rate were selected along the streamwise 
axis in the zone of established flow. Data planes slightly 
larger than the jet width were centered at these points and 
oriented perpendicular to -the streamwise axis . One hundred 
data points were selected in a symmetric square matrix with 
points most densely populated near the center. The planes 
were identified alphabetically and in consecutive order from 
"A" to "F", where "A" represented the plane nearest the 
nozzle. The positions were entered into the microcomputer 
and stored by plane on a floppy disk by the program LOAD_XYZ 
in Appendix B. Accompanying each data point was a probe 
deflection angle used to orient the probe parallel to the 
path of the jet to minimize interference. This angle was 
determined by evaluating the first derivatives of the 
equations developed for the streamwise axis and the half- 
widths and performing an interpolation based upon the data 
points' position relative to the two curves. 

C. DATA ACQUISITION 

The flow systems were placed into operation and the 
ambient and nozzle temperatures were monitored with the 
program T_SUBS to evaluate system stability and readiness 
for data acquisition. The system usually took approximately 
two hours to come into equilibrium. This could be monitored 
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by watching the jet nozzle temperature. When conditions 
were stable, data consisting of two hundred probe, ten 
ambient and ten nozzle temperature samples per position was 
collected, one plane at a time, by the program MAIN_T. 

The following information was stored on a floppy disk for 
each data point: x, y, and z coordinates; mean probe, 

nozzle and ambient temperatures and the standard deviation 
of the probe measurements. The data was transferred to 
the mainframe computer by using a modem, the microcomputer 
program SEND_DATA and the mainframe program GRAB. 

D. DATA REDUCTION 

The raw data was organized into a more usable format 
by the mainframe program TDATA which also converted the XYZ 
coordinates into the XSW system shown in Figure 2. The 
resulting data was selectively sent to the program C0NT0UR4 
which applied calibration coefficients to the temperature 
data and normalized it in the following manner; 




n a 



where T^ was the jet temperature as measured by the probe, 
T^ was the ambient fluid temperature and T^ was the 
temperature of the jet within the nozzle. 

Contour plots of this information, generated by the 
CONTOUR option of the graphics package DISSPLA [Ref. 8], 
are presented in Figures 30 through 35. 
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Figure 30 . Plane A Temperature Contour Plot 
(large scale) 
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Figure 31. Plane A Temperature Contour Plot 
(small scale) 
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Figure 32. Plane B Temperature Contour Plot 
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Figure 33 . 



Plane C Temperature Contour Plot 
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Figure 34 . 



Plane D Temperature Contour Plot 
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Figure 35. 



Plane E Temperature Contour Plot 



E. RESULTS 

The contour plots found in Figures 30 through 35 support 
the conjectures in Section II. C concerning the effects of 
the crossflow. Figure 30, a plot of plane A centered on 
the streamwise axis 46 degrees from horizontal and 7.327 mm 
(.288 in) downstream of the nozzle shows significant distor- 
tion. Figures 31 through 35 show planes A through E sequen- 
tially plotted on the same scale in order to observe overall 
jet behavior. Plane E was located 87.313 mm (3.438 in) 
downstream of the nozzle and 86 degrees from horizontal. 

It can be seen that as the jet traveled further downstream, 
the distorting effect grew progressively less, as expected. 

The rate of heat transfer from the jet to the ambient 
was calculated for each plane utilizing the temperature 
distribution matrix generated in the program CONTOUR4 in 
Appendix C and the following relationship: 



Q 



m n 

I I 

i=l j=l 



c A . . 
P 13 



U. . 

13 



T. . 
13 



where p was the relative density of the jet, A^^ was the area 

of each matrix segment, was the velocity in each segment, 

c was the specific heat of the jet and T. . was the tempera- 
p 13 

ture in each segment. Velocity was measured along the 
streamwise axis by a laser Doppler velocimeter (LDV) . It 
appeared to be nearly constant in the region of the jet 
observed. The mean and standard deviation was 44.875 mm/s 
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and 1.8 mm/s respectively. By assuming the Gaussian profile 
shown in Section II. B, velocity was determined at radii 
corresponding to each segment in the matrix mentioned 
above. The rates of heat transfer are shown in Table 2. 
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V. CONCLUSIONS AND RECOMMENDATIONS 



The objective of this thesis was to develop a computer- 
aided data acquisition system and construct a microthermo- 
couple probe to be used by follow-on students to study 
temperature distributions in turbulent buoyant jets. Sample 
data was taken to verify system operability. Based on 
results, the system performed in a satisfactory manner and 
will be an invaluable tool for subsequent studies of buoyant 
jets in a crossflowing ambient. 

Data point positions were hand calculated and loaded 
into the microcomputer with the program LOAD_XYZ. This was 
an extremely time consuming task and distracted the user 
from defining more than 100 data points per plane. The 
system can be greatly improved with the addition of micro- 
computer software that would automatically determine and 
load data point positions. Data point population could then 
be increased with ease which should result in smoother con- 
tour plots and more accurate heat transfer calculations. 

It was necessary to continuously add fresh water to the 
crossflow system as an equal amount was drained from it in 
order to maintain the crossflowing ambient at a constant 
temperature. This was because the cooling coil located 
within the flow settling chamber was inadequate to compen- 
sate for the heat added to the system by the jet and the 
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crossflow circulation pump by itself. Although it was 
possible to maintain the ambient constant within 1.4 C by 
this method, in the interest of conserving water, it is 
recommended that either a larger capacity chilled water 
bath or a cooling system that circulates a refrigerant 
rather than water be appropriated for this purpose. 
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APPENDIX A 



UNCERTAINTY ANALYSIS 

Experimental uncertainty was analyzed in accordance with 
the guidelines set forth by Holman [Ref. 9] . Uncertainties 
in the primary measurements, based upon manufacturer speci- 
fications and/or the number of significant digits which 
could be read, follow: 



1. 


Time (rotometer calibration) : 


.01 s 


2. 


Volume (rotometer calibration) : 


1.0 ml 


3. 


Rotometer reading: 


1.0 percent 


4. 


3-D positioning platform resolution: 


.1524 mm 


5. 


Probe deflection angle resolution: 


. 5 degrees 


6 . 


Thermocouple resolution: 


.0099 °C 


7. 


Thermocouple time constant: 


.004 s 


8. 


Time (crossflow velocity) : 


.01 s 


9 . 


Length (crossflow velocity) : 


.05 m 



Based on the above, the uncertainty of the flow rate of the 
jet was estimated to be 1.0 ml/s. Maximum uncertainty in 
the position of the probe's tip due to the uncertainty in 
deflection angle was determined to be .983 mm (.0387 in). 
Combined with the resolution of the positioning platform, 
the tip of the probe was positioned with an uncertainty of 
1.135 mm (.0447 in) in each plane, or with an overall 
uncertainty of 1.605 mm (.063 in) in three-dimensional 
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space. Temperature was measured at the approximate rate of 
100 samples per second, well within the constraints of the 
thermocouple time constant . Uncertainty in the temperature 
measurements were governed by the resolution of the analog- 
to-digital converter which was .0099 °C. Uncertainty in 
the crossflow velocity was .00625 m/s (.0205 ft/sec). 
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APPENDIX B 



MICROCOMPUTER PROGRAMS 

1. MAIN T 
10 

20 ! 

30 ? 

40 ! 

50 ! 

bO ! 

70 ! 

80 ! 

90 ! 

100 ! 

no ! 

120 ! 

130 r 

140 ! 

150 ! 

160 ! 

170 ! 

180 ! 

190 ! 

200 ! 

210 ! 

220 ! 

230 ! 

240 ? 

250 ? 

260 ! 

270 ! 

280 ! 

290 f 

300 Length _probe*4 , 4375 

310 Leng t h_arm*20 . 0 

320 ! 

330 OPTION BASE 1 

340 DIM Coef ( 12) ,X<500) . Y<500> ,Z(500) .P robe_ang le(500 ) 
350 LOADSUB ALL FROM “T_SUBS** 

360 LOADSUB ALL FROM *’MTR_SUBS” 

370 LOADSUB ALL FROM “PROBE.SUBS” 

380 CALL Retr ieve_coef (Coef < *) , ‘*mo tor_coef ’‘ ) 

390 ! 



400 ! 


2. Input desired positions from disk 


410 ! 






420 


I-l 




430 


INPUT 


"Angle of arm relative to + Y-ax is^‘* , Ang le_arm 
"Filename for positioning data'?" , F 1 1 ename'i 


440 


INPUT 


450 


ASSIGN 


‘?File4 TO Filenames 


460 Go on 




» 


470 




ENTER C«iFiie4;X(I) .Y(I) ,Z(I) 


480 ' 




1 


490 




ir X(I)O-100 THEN 


500 




IF X(I)=-999 THEN 


510 




Probe_ang le ( I ) = Y < I ) 


520 




P angie= Y (I) 


530 




ELSE 


540 




Probe angie(I)=P angle 


550 




1 = 1 + 1 


560 




END IF 


570 




f 


580 




GOTO Go on 


590 




END IF 


600 




» 



MAIN^T 

This program coordinates the entire 
data-taking evolation for measuring 
temperature distributions in buoyant 
jets . 



1 . Load all subprograms 

2. Input desired positions from a 
disk file of the form : "RUNEX** 
(value ’’XX" IS the run number) 

3. Move the 3-D positioning 
platform to each position 

4. Align the thermocouple probe un th 
the jet streamuise axis 

5. Obtain 200 temperatures at each 
position and compute the mean and 
standard deviation <sd) 

6 . Write to disk: 

a. T(C) 

b . sd 

c. x.y,i in (mm) 

nozrie centerline is (0,0.0) 

Probe and arm dimensions in inches: 

a. Shortest arm length * 20.0 

b. Longest arm length * 24.0 
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620 

630 

6^0 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

e_ang 

840 

850 

860 

870 

380 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 

1040 

1050 

1060 

1070 

1080 

1090 

1100 

1110 

1120 



GOTO Go_on 
END IF 

f 

ASSIGN C^File4 TO * 

N 1 terns* I- 1 

f 

!3. Begin loop to take data at each point 

t 

INPUT ••NAME OF FILE WHERE DATA IS TO BE ST0RED7*\Fi lenajne 1 $ 
Records* < N i teins*8*7/256 ) +2 
CREATE BOAT F i lenajne 1 2 , Records 
ASSIGN S»Filel TO Filena/nelS 



FOR I _position-l TO Nitems 

a. Move milling machine and move the 
thermocouple probe. 

I-I ^position 

CALL Move_ldv_to< X < I) ,Y(I) ,Z(I) ,Lengh t_arm »Length_probe , Angle_arm 
1 e( I) , Coef ( * ) ) 

I 

? 

!b. Obtain mean temperature and sd 

f 

CALL T_couple<T /’PROBE’*. “C*\ 200, St.dev) 

f 

!c* Measure ambient and jet temps 

f 

CALL T_coup le< T_amb lent , ’’AMBIENT” , ”C” , 1 0 , Sd) 

CALL T_coup le< T_noz z le , "NOZZLE” ,”C” , 1 0 ,Sd> 

I 

!d. Write ail information to disk 

I 

OUTPUT ^Filel ;X<I) ,Y(I) ,Z<I) 

OUTPUT ^F 1 1 e 1 ; T , T_amb i en t , T_no z z 1 e 
OUTPUT •^Filel :St_dev 
NEXT Imposition 

f 

?4, Close files 
? 

OUTPUT ^Filel ;-100 
ASSIGN @File1 TO * 

I 

!5. Term 1 nate ‘program 

f 

PRINT "All done!" 

BEEP 

BEEP 

END 
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2 . 



PROBE SUBS 



10 

20 

30 

AO 

50 

60 

70 

80 

90 

100 

no 

120 

130 

140 

150 

160 

170 

180 

ISO 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 



! PROBE^SIJBS 

I 

! This prograim moves the temperature probe 
! to desired angles of deflection. 

I 

! NOTE: Calibration coefficients are 

! entered in SUB Read_angle, 

! beginning at line 290. 

f 

CALL Read_angle< Angle) 

PRINT ’’The probe is presently at”;Angle;" degrees from horizontal. 
INPUT "What IS your desired angle for the probe'?" , Des i red_angle 
CALL Probe_move < Des 1 red_ang 1 e ) 

GOTO 120 
END 



SUB Read_angle< Actual_ang le) 

!This program reads the present angle 
!of the probe and returns it. 

! 0 degrees * horizontal 

!90 degrees = vertical » dounward 



AO » A1 and A2 “ coefficients for a 
second order curve fit of mV vs angle 
of deflection data. 



Coefficients for 8 June data follow: 



A0»-23. 8657824056449 
A1=3.09995506283164E-2 
A2=-l .22922305456149E-6 
! ? FORMAT A/D CARD 

OUTPUT 723:"CCJT" !Clear Relay Card 
OUTPUT 723;"CC»3T" IClear A/D Card 
OUTPUT 723;"CC»7T" IClear Digital Card 
OUTPUT 723;"SF,3,3.3,1 .25, 12T" 

OUTPUT 723;"0B J , 10 JT" ! CLOSE RELAY 
OUTPUT 723;"IP,3T" ! START A/D 

ENTER 72301 ;V 

Actual_angle=A0+A1*V+A2*V*V 

SUBEND 



SUB Probe_move < Des 1 red_ang 1 e ) 

tThis subprogram moves the probe to the 
Idesired angle 

! 0. Check to see if the angle is in 

! an acceptable range. 

f 

OUTPUT 723;"QP. I »0T" IClear Relays 
CALL Clear_screen 

I 

IF Desi red_ang 1 e>90 THEN 
BEEP 3400.1 
BEEP 1700,1 
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600 




BEEP 3400,1 


610 




PRINT "Desired angle exceeds 90 degrees 


620 




SUBEXIT 


630 




END IF 


6^0 




f 


650 




IF Desired angle<0 THEN 


660 




BEEP 3400,1 


670 




BEEP 3800,1 


680 




BEEP 3600.1 


690 




PRINT "Desired angle is negative???" 


700 




SUBEXIT 


710 




END IF 


720 


f 




730 


t 


1 .a. Clear the digital output card. 


7^0 


f 


b. Format the A/D card. 


750 


t 


c. Close the relay that corrects 


760 


f 


the probe potentiometer to the 


770 


f 


A/D converter. 


780 


t 




790 




OUTPUT 723;"OP.7,OT" 


800 




OUTPUT 723;"SF,3.3,3 J .25.3T" 


810 




OUTPUT 723;"CC,1T" 


820 




OUTPUT 723;"0B,1 ,10, IT" 


830 


t 




8^0 


\ 


2. Define the acceptable tolerance 


850 


\ 


in the angle (degrees) . 


860 


t 




870 




Tolerance- . 5 


880 


1 




890 


f 


3. Control loop. 


900 


t 




910 Repeat: 


CALL Read angleCActual angle) 


920 




PRINT "ANGLE 


930 




PRINT USING "DDD,DD":Actual_angle 


9^0 




BEEP Ac tual_ang le* 1 00 , . 05 


950 




Angle e r ror =• ( Des i red angle-Actual angle) 


960 




1 


970 




IF ABS(Angle er ror ) >To lerance THEN 


980 




t 


990 




IF Angle_er ro r >*0 THEN 


1000 




Direct ion$-"Doun" 


1010 




ELSE 


1020 




D 1 rec t ion$-"Up" 


1030 




END IF 


1040 




1 


1050 




CALL Motor go(Di rectionS) 


1060 




f 


1070 




GOTO Repeat 


1080 




END IF 


1090 


OUTPUT 


723;"OP,7,OT" 


1100 


OUTPUT 


723;"CC,1T" IClear relay card. 


I 110 


SUBEND 




1120 


j 




1130 






1140 


j 




1150 


SUB Motor_go<Di rect ion$) 


1160 


IF 


D i rec t ion$-"Up" THEN 


1 170 




Lbi t = 8 


1180 


END 


1 IF 


1190 


j 
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1200 

1210 

1220 



Lb 1 1 ^7 
END IF 



IF D 1 rec t ionS="Doun” THEN 



1230 ! 

1240 OUTPUT 723;**QP,7’*;2''Lbit;“T" 

1250 SUBEND 
1260 ! 

1270 ! 

1280 ! 

1290 SUB Clear_screen 
1300 I Clear the CRT. 

1310 ! 

1320 OUTPUT 2 USING ’V , B** ; 255 , 75 

1330 GCLEAR 

1340 SUBEND 
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3. MTR SUBS 



10 ! MTRJ5UBS 

20 ! _ 

30 ! The follouing series of suoroutines 

^0 ! are utilized to caliDrate tne positioning 

60 ! platform and ultimately to move the probe 

70 ! tip to desired positions uithin tne jet. 

80 ! 

SO SUB Drau^flume 

100 ! 

MO ! Draw the Buoyant Jet Flume on the CRT. 

120 ! 

130 GCLEAR 

140 GRAPHICS ON 

150 WINDOW 0.48,0,38 

160 LINE TYPE 1 

170 MOVE 6,6 

180 ! Drau the top viey 

ISO IDRAW 36.0 

200 IDRAW 0,10 

210 IDRAW -36,0 

220 IDRAW 0.-10 

230 ! 

240 IMQVE 6,0 

250 IDRAW 0,10 

260 IMOVE 24,0 

270 IDRAW 0,-10 

280 1 label ’*Top'\ 

290 MOVE 0.10 

300 CSIZE 5 

310 LABEL “Top” 

320 ! Drau the side vieu. 

330 MOVE 6,22 

340 IDRAW 36,0 

350 IDRAW 0,13 

360 IDRAW -36,0 

370 IDRAW 0,-13 

380 IMOVE 6,0 

390 IDRAW 0,13 

400 IMOVE 24,0 

410 IDRAW 0,-13 

420 ! label “Side” 

430 CSIZE 5 

440 MOVE 0.29 

450 LABEL “Side” 

460 ! 

470 ! Label the picture. 

480 ! 

490 MOVE 11,35 

500 CSIZE 7 

510 LABEL “BUOYANT JET FLUME” 

520 ! 

530 ! Put on the nozzle. 

540 ! 

550 MOVE 14.22 

560 IDRAW 0,2 

570 IDRAW .25,0 

530 IDRAW 0,-2 

590 ! 

600 IMOVE 0,-2 

610 CSIZE 3 
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620 


LABEL 


”no z z 1 e" 


630 


1 




6^0 




14.11 


650 


CSIZE 


3 


660 


LABEL 


”o’* 


670 


f 




680 


! Ind 


icate the direction of flow. 


690 


f 


* 


700 


fio^E 


B J 1 


710 


IDRAW 


2.0 


720 


IDRAW 


-.5. -.5 


730 


IDRAW 


.5, .5 


7^0 


IDRAW 


-.5. .5 


750 


t 




760 


fio^E 


<M 

CD 


770 


IDRAW 


2,0 


780 


IDRAW 


-.5. -.5 


790 


IDRAW 


.5. .5 


800 


IDRAW 


-.5. .5 


810 


f 




820 


SUBENO 


830 


f 




8^0 


» 




850 


» 




860 


» 




870 


f 




880 


f 




890 


SUB C 


alibrate<Fi lenajneS) 


900 


1 




910 


OPTION BASE 1 


920 


DIM 


Coef (12) 


930 


f 




9A0 


! Calibrate the positioners on the 


950 


! m 


illing machine movement. 


960 


! Q 


nto disk, urite out the calibration 


970 


! coefficients and the hard boundaries 


980 


! that must be observed! 


990 


! This file uill be called ’‘Motor coef”. 


1000 


1 




1010 


! A. 


Position the probe volume at the 


1020 


1 


wall of the tip of the nozzle. This 


1030 


f 


position IS (0,0.0). All readings 


10^0 


t 


uill be in inches. Read all three 


1050 


f 


potentiometers. Ask the user for the 


1060 


1 


nozzle outer diameter and compute the 


1070 


\ 


zero position. Ask the user for the 


1080 


1 


milling machine readings. 


1090 


t 


1100 


! B. 


Next, move the bed to some new posi- 


1 1 10 


» 


tion using the override switches. 


1120 


» 


Take readings from the pots and ask 


1130 


f 


for the milling machine readings. 


11^0 


I 


Compute the calibration coefficients. 


1150 


t 




1 160 


! C. 


Move the bed to each of the extremes 


1 1 70 


f 


in the X, Y. and Z directions using 


1 180 


I 


the override switches and have the 


1190 


f 


user tell the computer when each of 


1200 


i 


these boundaries are hit. Enter each 


1210 


1 


of these onto the disk file. 
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1220 

1230 

1240 

1250 

1260 

1270 

1280 

1290 

1300 

1310 

1320 

1330 

1340 

1350 

1360 

1370 

1380 

1390 

1400 

1410 

1420 

1430 

1440 

1450 

1460 

1470 

1480 

1490 

1500 

1510 

1520 

1530 

1540 

1550 

1560 

1570 

1580 

1590 

1600 

1610 

1620 

1630 

1640 

1650 

1660 

1670 

1680 

1690 

1700 

1710 

1720 

1730 

1 740 

1750 

1760 

1770 

1780 

1790 

1800 

1810 



D. Disk file ”Mo tor_coef : 

1 . x_zero , x_s lope 

2. y_zero, y_slope 

3. z_zero, z_slope 

4 . x_m 1 n , x_max 

5. y_nu n ♦ /_max 

6. z_max 

BEEP 

PRINTER IS 1 
GCLEAR 

OUTPUT 2 USING ’V . B’* ; 255 . 75 

f 

PRINT •*! 

PRINT •* 

PRINT 
PRINT •* 

PRINT •• 

PRINT ” 

PRINT •• 

PRINT ” 

PRINT 
PRINT ” 

PRINT 
PRINT '• 

PRINT 
PRINT ’• 

PRINT •* 

PRINT 
PAUSE 
BEEP 1500 , .1 

INPUT “1. Nozzle Q.D. ( i nches ) ?’* , No z z le_od 
BEEP 2000 . . 1 

INPUT ”2. X (in), {+ into flame)?’\X_l 
BEEP 2500, .1 

INPUT ”3. y (in), (+ along flume to the right}?’*, Y_1 
BEEP 3000 , . 1 

INPUT -4, z (in), upyard}?”,Z_l 

I 

(:ALL Read^ot(”X”,Vx_l ) 

CALL Read_j)ot(”Y”.Vy_l ) 

CALL Read _pot <”Z” , Vz^l ) 

f 

I 

CALL Clear_screen 

PRINT ”C. Move the milling machine to a neu” 

PRINT ” position in 3-D, by at least 5 
PRINT ” inches in each direction. ” 

BEEP 3200. . 1 

INPUT -X. Y, Z in inches?” , X_2 , Y_2 , Z^2 
CALL Read^otC’X” .Vx_2) 

CALL Read _po t < ”Y” , Vy_2 ) 

CALL Read _po t < ”Z” . U z^2 ) 

I 

! C. Calculate the calibration coefficients 

f 

X_zero--( (Nozzle_od/2)+Vx_l*( (X_2-X_l ) / (Ux_2-Vx_l ) ) ) 
X_slope-<X_.2-X_l )/(Ux_2-Vx_l ) 



. POTENTIOMETER CALIBRATION: ” 

NOTE: 1. Probe must be horizontal” 

2. Arm must be parallel to 
the bed axis 

3. ALPHA - BETA ” 

A. Using the override switches, 
position the probe volume at 

the outer ual 1 of the tip of ” 
the nozzle. 

B. I will need the nozzle Q.D. and” 
the milling machine positon. ” 

M 

Hit <cont>> ” 
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1820 

1830 

18^0 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

1940 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

2040 

2050 

2060 

2070 

2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 

2300 

2310 

2320 

2330 

2340 

2350 

2360 

2370 

2380 

2390 

2400 

2410 



I 



Y_zero»-My_l^(Y_2-Y_1 )/(Vy_2-Vy_l ) 
Y_slope-( Y_2-Y_l )/(Vy_2-Vy_l ) 

♦ 

Lzero--^z_l*(Z_2-Z_l )/(Vz_2-Vz_1 ) 
Z_slope-(Z_2-Z_l )/(Vz_2-Vz_l ) 



machine to the 
' X ' . (Towards 
<con t> . 



machine to the 
' y ' . Hit <cont> 
the flume} 



D. Find the physical boundaries for each 
di rec t ion . 

CALL Clear_screen 

PRINT ‘*1. Move the milling machine to the 
PRINT ** minimum value of 'x'. Hit <cont> 
PRINT '*(Away from the flume, backwards/ 
PAUSE 

CALL Read_pot < "X** , V) 

X_m in = X_ze ro + X_s lope**^ 

PRINT "2. Move the milling 
PRINT “ maximum value of 
PRINT •• the flume}. Hit 
PAUSE 

CALL Read_po t ( “X" , U ) 

X_max = X_ze ro + X_s lope**^ 

PRINT "3. Move the milling 
PRINT " minimum value of 
PRINT "(To the left along 
PAUSE 

CALL Read_pot ("Y" , U) 

Y_min= Y_zero+ Y_s lope*V 
PRINT "4. Move the milling machine to the 
PRINT maximum value of 'y'. Hit <cont> 

PAUSE 

CALL Read ^ot(”Y**.V) 

Y_max =- Y_zero + Y_slope«U 
PRINT '*5. Move the milling 
PRINT " minimum value of 
PRINT "(Downwards}'* 

PAUSE 

CALL Read_j>o t ( "Z** , V) 

Z_m i n = Z_zero + Z_s lope*V 
PRINT "6. Move the milling 
PRINT '* maximum value of 
PAUSE 

CALL Read^ot(**Z**,V) 

Z_max*Z_zero+Z_slope*U 

I 

! E. Write out the file "Motor_ 

j 

ON ERROR GOTO Purge^file 
Reenter: CREATE BOAT Filenames,! 

ASSIGN '^File TO Filenames 



machine to the 
' z ' . Hit <con t > 



machine to the 
' z ' . Hit <con t > 



coef * 



OUTPUT ?lFile;X_zero.X_slope 
OUTPUT »^F 1 le ; Y_zero . Y_s lope 
OUTPUT ^File:Z_zero,Z^slope 



OUTPUT '?F 1 le ; X_min , X_max 
OUTPUT SF 1 1 e : Y_m I n , Y_max 
OUTPUT ?iF 1 1 e ; Z_m 1 n , Z_max 
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2^20 ! 

2430 ASSIGN ?File TO • 

2440 SUBEXIT 

2450 Purge_file: PURGE Filenames 

2460 GOTO Reenter 

2470 ! 

2480 SUBEND 
2490 ! 

2500 ! 

2510 ! 

2520 ! 

2530 ! 

2540 ! 

2550 SUB Read_pot<DirectionS, Value) 

2560 ! 

2570 ? 

2580 ! Read one potentiometer and return a volt- 

2590 ! age. 

2600 ! 

2610 ! R3 Pot X 

2620 ! R4 -- Pot Y 

2630 ! R5 -- Pot Z 

2640 ! 

2650 ! 

2660 IF Directin$ = **X** THEN Relay-3 
2670 IF D i rect ion$-‘'Y*’ THEN Relay*4 
2680 IF Direction$'”Z‘* THEN Relay-5 
2690 ! 

2700 ! 

2710 OUTPUT 723:"0P J .OT** 

2720 OUTPUT 723 ; ”0B J , " ; Re lay ; ’M T‘* 

2730 ! 

2740 OUTPUT 723;”IP,3T** 

2750 ENTER 72301 ;Value 
2760 ! 

2770 ! 

2780 SUBEND 
2790 ! 

2800 ! 

2810 ! 

2820 ! 

2830 ! 

2840 ! 

2850 SUB Clear_screen 
2860 ! Clear the CRT, 

2870 ! 

2880 OUTPUT 2 USING ‘V , B’* ;255 . 75 

2890 GCLEAR 

2900 SUBEND 
2910 ! 

2920 ! 

2930 ! 

2940 ! 

2950 ! 

2960 ! 

2970 SUB Mo tor ( D 1 rec t lonS , Ro tat lonS) 

2980 ! 

2990 ! Turn on the motor in the requested direc- 

3000 ! tion (x,y,z) yith the requestion rota- 

3010 ! tion (CU. CCW) . 
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3020 ! 

3030 D i rS=D 1 rec t lonS 

30^0 Ro tS-Rotat ion$ 

3050 IF DirS = ”X’* AND Rot$-’*CU” THEN Lbit-2 

3060 IF DirS-”X” AND RotS-‘*CCW‘ THEN Lbit-1 

3070 IF DirS-”Y” AND RotS-’*CU” THEN Lbit = 4 

3080 IF Dir$ = '*Y’* AND Rot$ = ”CCW'* THEN Lbit»3 

3090 IF DirS=”Z” AND Rot$**TCW" THEN Lbit-5 

3100 IF Dir$-”Z‘* AND Rot$ = ’*CW’ THEN Lbit-6 
3110 I 

3120 I 

3 1 30 OUTPUT 723 ; •*0P , 7 ; 2"Lb 1 1 ; ’’T” 

31^0 SUBEND 
3150 ! 

3160 ! 

3170 ! 

3180 ! 

3190 ! 

3200 ! 

3210 SUB Motor_stop 
3220 ! 

3230 ! Stop ail motors! 

32^0 ! 

3250 OUTPUT 723 :‘’0P . 7 ,or* 

3260 SUBEND 
3270 ! 

3280 ! 

3290 ! 

3300 ! 

3310 ! 

3320 ! 

3330 SUB Retr ieve_coef (Coef ( *) ,Fi lena/neS) 

33^0 ! 

3350 OPTION BASE 1 

3360 ! 

3370 ! Retrieve the potentiometer calibration 

3380 ! coefficients from a disk file called 

3390 ! "No tor_coef . Place these in an array. 

3^00 ! 

3^10 ASSIGN iiFile TO Filenames 

3^20 ! 

3^30 FOR I-l TO 12 STEP 2 
3440 ENTER «File;Coef(I) »Coef<I+1) 

3450 NEXT I 

3460 ! 

3470 ASSIGN @File TO * 

3480 SUBEND 
3490 ! 

3500 ! 

3510 ! 

3520 ! 

3530 ! 

3540 ! 

3550 SUB Label^oint<X,Y,Z, Symbols) 

3560 ! 

3570 ! Label a point on the Flume diagram 

3580 ! using the symbol specified. 

3590 ! 

3600 PRINT TABXY(1J3) ."(X.Y.Z) inches - 
3610 PRINT USING "DDD . ODD" ; X . Y , Z 
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3620 

3630 

3640 

3650 

3660 

3670 

3680 

3690 

3700 

3710 

3720 

3730 

3740 

3750 

3760 

3770 

3780 

3790 

3800 

3810 

3820 

3830 

3840 

3850 

3860 

3870 

3880 

3890 

3900 

3910 

3920 

3930 

3940 

3950 

3960 

3970 

3980 

3990 

4000 

4010 

4020 

) ) 

4030 

4040 

4050 

4060 

4070 

4080 

4090 

4100 

4110 

4120 

4130 

4140 

4150 

4160 

4170 

4180 

4190 

4200 



ir SymbolS-’*X” THEN 

MOVE 14+Y.24+Z 
IMQVE -.2. -.2 
IDRAW 0. .4 
IDRAW .4,0 
IDRAW 0,-.4 
IDRAW -.4,0 

I 

MOVE 14+Y,11>X 
IMQVE -.2, -.2 
IDRAW 0, .4 
IDRAW .4,0 
IDRAW 0,-.4 
IDRAW -.4.0 
END IF 
\ 

IF Symbols-”-*-” THEN 

j 

MOVE 14 + Y, 1 1^-X 
IMOVE 0,.2 
IDRAW -.2, -.4 
IDRAW .4,0 
IDRAW -.2, .4 

I 

MOVE 14+Y.24+Z 
IMOVE 0, .2 
IDRAW -.2, -.4 
IDRAW .4,0 
IDRAW -.2. .4 
END IF 
SU8END 



SUB Move__ldv_to( X , Y , Z , Leng th_arm , Leng th__p robe , Ang I e_arm .Angle _p robe , 
DEG 

j 

OUTPUT 723;”CC,1T” 

OUTPUT 723;”SF,3,3.3.1 .25.12T” 

PRINTER IS 1 

ON KEY 0 LABEL ”ABORT” CALL Stop_alI 
OPTION BASE 1 
Tolerance* . 006 

I 

! Move the probe to the position indicated 
I in (inches) relative to the noizle tip. 

» 

! Dimensions are in inches and degrees 
» 

! A.1 Load in the calibration coefficients. 

I 

X_z e ro-Coef ( 1 ) 

X_s lope = Coef ( 2 > 
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^210 Y_:ero=Coef (3 ) 

4220 Y_slope-Coef (4) 

4230 Z_zero=Coef ( 5 ) 

4240 Z_s lope=Coef (6 ) 

4250 X_m i n*Coef ( 7 ) 

4260 X_max =Coef (8 > 

4270 Y_m i n *Coef ( 9 > 

4280 Y_max*Coef n 0 ) 

4290 Z_m i n=Coef ( 11 ) 

4300 Z_max -Coef < 1 2 ) 

4310 ! 

4320 ! 8.2 Move the probe 

4330 ! 

4340 CALL Probe_moe( Ang le_probe ) 

4350 ! 

4360 ! B, Does (X,Y,Z) lie uithin the per- 

4370 1 mitted boundaries'? 

4380 ! 

4390 IF X>X_max OR X<X_rrun OR Y>Y_max OR Y<Y_min OR Z>Z_max OR Z<Z_nin THEN 



4400 

4410 

4420 

4430 

4440 

4450 

4460 

4470 

4480 

4490 

4500 

4510 

4520 

4530 

4540 

4550 

4560 

4570 

4580 

4590 

4600 

4610 

4620 

4630 

4640 

4650 

4660 

Ang le_p 

4670 

Ang le_p 

4680 

Ang le_p 

4690 

4700 

4710 

4720 

4730 

4740 

4750 

4760 



BEEP 1 700 , .5 
BEEP 2000, .5 
BEEP 1 700, .5 
BEEP 2000, .5 

PRINT '’Desired point is out of range!'* 

SUBEXIT 
END IF 

C. Find out where the probe is now, draw 
the flume on the CRT, and label the 
desired position. 

1. Sound warning, movement immenent! 

CALL Clear_screen 

PRINT TABXY( 1 , 12) .“MOVEMENT OF MILLING MACHINE IMMINENT! !!“ 

FOR 1=1 TO 4 
BEEP 1200, .1 
BEEP 1 700, .1 
BEEP 2200,. 1 
BEEP 2700, .1 
NEXT I 

CALL Clear_screen 

\ 

\ 

OUTPUT 723;“CC,1T“ ! CLEAR RELAY CARD. 

CALL Posi t ion( “X” ,X_actual , Values haf t ,Length_arm , Leng th_probe .Ang le_arm . 
robe , Coef ( *) ) 

CALL Pos 1 1 ion< ’*Y“ , Y_actual , Values haf t , Length_arm ,Leng th_probe . Ang le_arn . 
robe , Coef ( *) ) 

CALL Pos 1 1 ion< “Z“', Z_actual , Val u_shaf t .Leng th_arm , Leng th_p robe . Ang le_arn . 
robe , Coef ( *) ) 

t 

f 

CALL Draw_flume 

CALL Label_po int(X_actual , Y_actual , Z_actual , 

CALL Label_point(X.Y,Z,'’X") 

D. Move each motor to bring the error 
between actual and desired position 
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4770 

4780 

4790 

4800 

4810 

4820 

4830 

4840 

4850 



into tolerance. 



X_olcl*X_actLiai 
X_node : Xe r ror =X-X_ac t ual 

IF ABS( Xerror ) >Tolerance THEN 
IF Xerror>0 THEN Rot$-‘TCW” 

IF Xerror<0 THEN Rot$"**CU’* 

CALL Motor ( ”X" , Rots ) 

CALL Posi t ion< “X'* , X_ac taal , Valu_shaf t , Leng th_arm , Length_probe , Ang I e_ai 



Ang le_probe , Coef ( • ) ) 



4860 

al ) 

4870 

4880 

4890 

4900 

4910 

4920 

4930 

4940 

4950 

4960 

4970 

4980 



CALL P lo t_path (X_o Id , Y_ac tual , Z_ac tual , X_ac t ual , Y_actual , Z_a< 



X_o ld’X_ac taal 
GOTO X_node 
END IF 



Y node 



Y_o ld“Y_actaal 
Y e r r o r ’ Y - Y _ ac t aa 1 
IF ABS( Yerror )>Tolerance THEN 
IF Yerror>0 THEN Rot$= 

IF Yerror<0 THEN Rot$= 

CALL MotorC*Y”,RotS> 

CALL Pos 1 1 ion< "Y** , Y_actual , Valu_shaf t ,Length_arm » Length_p robe , Ang le_a 



*CCW 

•cw 



Angle _probe , Coef ( • ) ) 



4990 

al) 

5000 

5010 

5020 

5030 

5040 

5050 

5060 

5070 

5080 

5090 

5100 

5110 



CALL P lo t_path< X_ac tual ♦ Y_o Id , Z_ac toal , X_ac taal , Y_ac toal , Z_a( 



Y_old* Y_actaal 
GOTO Y_node 
END IF 



Z_o ld*Z_actoal 
Z_node : Zerror’Z-Z_actoal 

IF ABS( Zerror)>Tolerance THEN 
IF Zerror>0 THEN RotS»'*CW’ 

IF Zerror<0 THEN RotS--CCW 
CALL Motor(“Z”.RotS) 

CALL Posit ion< **Z” , Z_actaal , Val a_shaf t , Length_arm ,Length_probe . Angle_a 



Ang le_p robe , Coef ( • > ) 



5120 

al ) 

5130 

5140 

5150 

5160 

5170 

5180 

5190 

5200 

5210 

5220 

5230 

5240 

5250 



CALL P lo t_path( X_ac toal , Y_actoal , Z_o Id , X_actaal , Y_actoal , Z_a 



Z_old»Z_actoal 
GOTO Z_node 
END IF 



CALL Motor_stop 



FOR I»1 TO 4 
BEEP 2400, .2 
BEEP 4800. .2 
NEXT I 



5260 SBEND 
5270 ? 



5280 

5290 

5300 
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5310 ! 

5320 ! 

5330 SUB Pos ition<DirectionS,Ualue ,Ualu_shaf t »Leng th_arm , Leng th_p robe . Ang le_arm . 
Ang le_probe , Coef ( * ) ) 

5340 ! 

5350 OPTION BASE I 
5360 DEG 
5370’ ! 

5380 ! ’ Return the position (inches) for the 

5390 ! appropriate direction relative to the 

5400 ! nozzle tip. 

5410 ! 

5420 CALL Read _po t ( D i rect lonS , Uo 1 tage) 

5430 ! 

5440 X_zero=Coef ( 1 ) 

5450 X_s 1 ope=Coef (2 ) 

5460 Y_zero=Coef <3 ) 

5470 Y_s lope-Coef ( 4 ) 

5480 Z_zero=Coef ( 5 ) 

5490 Z_s lope=“Coef (6) 

5500 ! 

5510 ! 

5520 D 1 rS»D i rec t ionS 

5530 ! 

5540 IF DirS-’*X*’ THEN 

5550 ^)al u_shaf t’X_zero + X_s lope*Uo 1 tage 

5560 Val ue*Ual u_shaf t-Leng th_arm*C0S ( Ang le_arm) »TAN( 45-Ang le_arm/2 . 0 ) 

5570 END IF 

5580 ! 

5590 IF Dir$-”Y’* THEN 

5600 Val u_shaf t’* Y_zero'»' Y_s lope*^)o I tage 

5610 Ual ue“Ual u_shaf t ■♦•Leng th _probe*< 1 . -CQS( Ang le _probe) ) +Length_arm^COS( 

Ang le_arm) 



5620 END IF 

5630 ! 

5640 IF DirS-**Z" THEN 

5650 Ual u_shaf t = Z_zero + Z_slope**Uol tage 

5660 Ual ue*Ualu_shaf t“Length_j5robe*SIN( Angle _probe) 

5670 END IF 

5680 ! 

5690 SUBEND 
5700 ? 

5710 ! 

5720 ! 

5730 I 
5740 ! 

5750 ! 

5760 SUB Stop_all 

5770 ! STOP ALL MOTORS AND QUIT 

5780 ! 

5790 GCLEAR 

5800 CALL Motor_stop 

5810 PRINT ’’MOTOR CONTROL ABORTED!!!” 

5820 PRINT ’’(HIT <C0NT> TO CONTINUE)” 

5830 PAUSE 

5840 SUBEND 
5850 ! 

5860 ! 

5870 ! 

5880 ! 
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5890 

5900 

5910 

5920 

5930 

5940 

5950 

5960 

5970 

5980 

5990 

6000 

6010 

6020 

6030 

6040 

6050 

6060 



SUB Plot ^ath<XKYl ,Z1 ,X2,Y2»Z2) 

f 

! Plot the path of the probe on the flame 
! diagram as the motors move the bed. 

! Lower plot followed by upper plot. 

» 

PRINT TABXYd J3) /*(X,Y,Z) inches - ; 

PRINT USING •’DD.DDD";X2,Y2,Z2 
MOVE 14+Yl »24+Z1 
DRAW 14+Y2.24+Z2 

f 

lioUE 14 + Yl ,1 M-Xl 
DRAW 14 + Y2.1UX2 

I 

SUBEND 
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4. T SUBS 



20 

30 

40 

50 

60 

70 

80 

90 

100 

110 

120 

130 

140 

150 

160 

170 

130 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

380 

390 

400 

410 

420 



460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

580 

590 

600 



! T.SUBS 

I 

I 

SUB T_coulp le ( Temperature , Cho iceS , ScaieS , No_read i ngs , S tdev ) 

SUBPROGRAM T_C0UPLE 

by Bill Cuibreth 
19 April 1984 

PURPOSE: This program is designed to 

read type T or E thermocouples and 
return the actual temperature to 
the calling routine. 

Temperature -- Temperature from the 
thermocouple in degrees F or C from 
the thermocouple identified by 
"Cho iceS” . 

Choices -- Thermocouple choice, current- 
ly: "AMBIENT", "NOZZLE", or. "PROBE". 

Scales -- "F" for Fahrenheit or "C" For 
Celsius, or "H" for histogram in "C" . 

No_readings -- Hou many readings of the 
same thermocouple should the routine 
taic e? 

Stdev -- The standard deviation of the 
temperature for the indicated number 
of readings in the units given by 
"Scales" . 



1. Open all relays and initialize the 
A/D converter. 

OUTPUT 723;"CC,1T" 

OUTPUT 723;"SF.3,3,3,1 .25,12T" 

t 

OUTPUT 723; "OP, 1 ,0T" 



2. Close the chosen relays. 



430 


f 


a. 


Ambient 


T — 


Type 


T. 


relays 


6,8 


440 


i 


b. 


Nozzle 


T — 


Type 


T, 


relays 


7,8 


450 


t 


c . 


ProDe 


T -- 


Type 


E, 


relay 


9. 



IF Scale$*"H" THEN 

! Plot a histogram using Celsius. 
H is togram$=" YES" 

Seal eS="C" 

ELSE 

H IS tog ram$“"N0" 

END IF 

» 

» 

IF ChoiceS“"AMBIENT" THEN 
TypeS="T" 

OUTPUT 723;"QB, 1 .6,1,8. IT" 

END IF 
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610 

620 

630 

640 

650 

660 

670 

680 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

340 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 



IF Choice$*’‘NGZZLE’* THEN 
TypeS = '*T'‘ 

OUTPUT 723;‘’0BJ ,7,1 ,3, ! T*‘ 

END IF 

t 

IF ChoiceS-’*PROBE” THEN 
TypeS^ "E" 

OUTPUT 723;”QB J ,9.1T** 

END IF 

» 

3, Taice an A/D conversion and convert 
into temperature. 

$um = 0 
Su/n I =0 

WAIT I 

IF HistogramS-'*YES” THEN GOSUB Set_up_histo 

t 

f 

FOR 1*1 TO No_readinqs 
OUTPUT 723;’*IP,3T’* 

ENTER 72301 ;A 

IF Type$-”T” THEN 
A*A/1000 
END IF 

IF TypeS-”E’* THEN 
A*A/1000 
END IF 

! PRINT ”V(mV) - *’;A 

GOSUB Convert_t 
! BEEP A*100..01 

Su/n*Su/n + A 
Sumi -Su/nl +A»A 

f 

IF H IS togram$*‘*YES’* THEN GOSUB Plot _point 



990 




NEXT I 


1000 


t 




1010 




Temperature“Sum/No_readings 


1020 




t 


1030 




f 


1040 




IF No_read 1 ngs- 1 THEN 


1050 




S tdev=*0 


1060 




ELSE 


1070 




$ tdev-SQR ( ABS ( < Sum 1 -No n 


1080 




END IF 


1090 


1 




1100 


i 


4. Open ail relays. 


1 110 


t 




1120 




OUTPUT 723;’*0P,1 ,0T‘* 


1130 


1 




1140 


SUBEXIT 


1150 


f 




1160 


» 




1170 


» 




1180 


Convert t: ? 


1190 




! This subroutine c 


1200 




! from a thermocoup 
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1210 
1220 
1230 
12^0 
1250 
1260 
1270 
1280 
1290 
1300 
1310 
1320 
1330 
13^0 
1350 
1360 
1370 
1380 
1390 
MOO 
MIO 
M20 
M30 
1 ^40 
1450 
1460 
1470 
1480 
1 490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1 700 
1710 
1 720 
1730 
1740 
1750 
1760 
1 770 
1780 
1790 
1800 



! erature. 



IF TypeS='T* THEN 

A-2.5661297E+MA-6. 1954869E-1*A*A^^2.2181644E-2*A^3-3.55009E-4*AM 

END IF 



IF TypeS-”E” THEN 

A-1 .7022525E + 1*A-2.2097240E-1*A*A+5.4809314E-3*A''3-‘5.7669892E-5*AM 
END IF 

t 

! Fix the scale. 

f 

IF ScaieS»"F” THEN A-1.8*A+32 
RETURN 



Set_LLp_histo : ! 

f 

! Set op a histogram of temperatare 
! versus number of counts. 

t 

! 1. Zero out the Height<*) array. 

f 

Dili Height<203) 

f 

FOR I-l TO 202 
He ight ( I ) -0 
NEXT I 



DUMP DEVICE IS 701 
GINIT 

OUTPUT 2 USING 'V , B” ; 255 » 75 

GRAPHICS ON 

FRAME 

WINDOW -100 JOO.-IO JOO 
MOVE -65,92 
CSIZE 7 

LABEL "TEMPERATURE HISTOGRAM" 

AXES 25,10,0,0,4,5,3 

PEN -1 

MOVE 0,-10 

DRAW 0,0 

MOVE 0,90 

DRAW 0,100 

PEN 1 

f 

! Take 10 temperature readings to get 
! the scale. 

f 

Sum*0 
Sum I =0 

FOR Ml TO 10 

OUTPUT 723;"IP,3T" 



83 



1810 




ENTER 72301 


1320 




NEXT I 


1830 




t 


18^0 




FOR I-l TO 100 


1850 




OUTPUT 723;**IP.3T“ 


1860 




ENTER 72301 ;A 


1870 




IF TypeS="T” THEN 


1880 




A*A/1000 


1890 




END IF 


1900 




IF Type5***E” THEN 


1910 




8*8/1000 


1920 




END IF 


1930 




GOSUB Convert_t 


1940 




Sujn*Siun + 8 


1950 




Sum 1 *Sum 1 +8»A 


1960 




NEXT I 


1970 




T mean*Sum/100 


1980 




Sd*SQR(8BS(Sum1-100^T mean'2)/99) 


1990 




f 


2000 




? Change the uindoui to extent from 


2010 




f -4*Sd to +4*Sd. 


2020 




t 


2030 




MOVE -12.-7 


2040 




CSIZE 4 


2050 




LABEL USING ’*DDD.DD**;T mean 


2060 




MOVE 15,-7 


2070 




LABEL ’’C” 


2080 




t 


2090 




MOVE -87,-7 


2100 




LABEL USING **DDD.DD”;T mean-3*Sd 


2110 




MOVE 63,-7 


2120 




LABEL USING ’*DDD.DD’*;T mean+3*Sd 


2130 




f 


2140 




1 


2150 




MOVE 5.47 


2160 




LABEL '’50** 


2170 




» 


2180 




f 


2190 




! 3. Calculate the window temper- 


2200 




? ature interval . 


2210 




t 


2220 




Interval *4*Sd/ 100 


2230 




f 


2240 




T_min*T_mean-4*Sd 


2250 




T_max=T mean+4»Sd 


2260 




1 


2270 




WINDOW T min.T max, -10. 100 


2280 


RETURN 


2290 


f 




2300 


f 




2310 


f 




2320 


f 




2330 


f 




2340 


t 




2350 


Plot point: ! 


2360 




» 


2370 




? Plot each temperature point as 


2380 




? received on the histogram. 


2390 






2400 




Freq*( (A-T_min)/( T_max - T_m in) )*3000+l 
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2^1 0 BEEP Freq ,.01 

2^20 J-0 

2^30 » 

2^40 FOR T-T_fnin TO T_max STEP Interval 

2A50 J-J+1 

2460 IF P<-T THEN 

2470 He iqh t < J ) »He igh t < J ) + 1 

2480 GOTO Continue 

2490 END IF 

2500 NEXT T 

2510 ! 

2520 ! Drau the point. 

2530 ? 

2540 Continue: ! 

2550 ? 

2560 ! 

2570 MOVE A.Height(J) 

2580 CSIZE 4 

2590 LABEL 

2600 ! 

2610 ! Label the number of points and 

2620 ! the temperature, 

2630 ! 

2640 PRINT TABXY(2,3) ;”T(C)-“;A 

2650 PRINT TABXY(2,5) ;**Sample ^”;I 

2660 ! 

2670 RETURN 
2680 ! 

2690 ? 

2700 ! 

2710 ? 

2720 ! 

2730 ! 

2740 SUBEND 
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5. 



T CAL 



! 0 ! T_CAL 

20 ! 

30 IThis program records the output from the 
40 lambient, probe and nozzle themocouples 
50 !for calibration purposes. The data is 

60 Iprinted as yell as transmitted to a data 

70 !file. 

80 ! 

30 IData is recorded in the following order: 

100 ! a. Ambient temperature 

110 ! b. Probe temperature 

120 ! c. Nozzle temperature 

130 ! 

140 IPlace the ambient, probe and nozzle 

150 ! thermocoup les into a bucket of warm water, 

160 lexecute this program, then add ice. 

170 ! 

180 ! 

190 LOADSUB ALL FROM •’T^SUBS’* 

200 ! 

210 ! Identify the BOAT file 

220 ! 

230 INPUT ’*NAME OF FILE WHERE DATA IS TO BE STORED'^** , Fi lename I S 
240 ! 

250 ! Identify the number of data points desired 

260 ! 

270 INPUT "NUMBER OF DATA POINTS DESIRED':^" , N i terns 
280 ! 

290 ! Identify the number of samples per 

300 ! thermocouple per data point 

310 ! 

320 INPUT "NUMBER OF SAMPLES PER DATA POIN J?" .Samp le 
330 ! 

340 ! Initialize a counter 

350 ! 

360 1=1 

370 ! 

380 ! Create a data file 

390 I 

400 Reco rds* < N 1 tems*8*3/256 > +2 
410 CREATE BDAT F i lenamel S , Records 
420 ASSIGN ^Filel TO FilenamelS 
430 ! 

440 ! Take data and print results 

450 ! 

460 CALL T_coup le< T_amb 1 en t , " AMB lENT" , "C" , Samp le , Sd > 

470 BEEP T*50. .05 

480 PRINT "TAMB,SD=":T_ambient ,Sd 

490 CALL T_coup I e < T , "PROBE" , "C" , Samp le , S t_dev ) 

500 BEEP T*50, .05 

510 PRINT "TPR0BE,SD=";T.St^dev 

520 CALL T_coup 1 e ( T_no z z 1 e , "NOZZLE" , "C" , Samp 1 e , Sd) 

530 BEEP T*50, .05 

540 PRINT "TNOZL.SD=";T_nozzle,Sd 

550 ? 

560 ! Send data to the data file 

570 ! 

580 OUTPUT ?File!;T ambient,T,T nozzle 
590 ! 

GOO ! Test to see if finished 
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bio 

620 

530 

6^0 

550 

650 

670 

680 

530 

700 

710 

720 

730 

7^0 

750 

760 

770 

780 

790 



IF I^Nitens THEN 
GOTO 750 
ELSE 

GOTO ^60 
END IF 

! Close the data file 

OUTPUT iFilel ;-100 
ASSIGN ^Fiiel TO » 

! Alert the user... the job is completed 

PRINT “ALL DONE*' 

BEEP 

BEEP 

END 



PROBE CAL 



6 . 



10 

20 

30 

^0 

50 

60 

70 

80 

90 

100 

no 

120 

130 

1^0 

150 

160 

170 

180 

190 

200 

210 

220 

230 

240 

250 

260 

270 

280 

290 

300 

310 

320 

330 

340 

350 

360 

370 

330 

390 

400 

410 

420 

430 

440 

450 

460 

470 

480 

490 

500 

510 

520 

530 

540 

550 

560 

570 

530 

590 

600 



probej:al 

This program uill provide the data 
necessary to calibrate probe deflection 
if used in the following manner: 

a. Switch OFF probe actuation power 
at the probe base. 

b. Swing the probe arm out of the 
tank such that the probe is 
accessab le , 

c. Attach the calibration panel to the 
probe assembly, taking care not to 
damage the glass probe. 

d. Run the program.. .it will ask for 
the desired position in mV. The 
following guidelines apply: 

I. If it IS desired to lower 
the probe, type the extreme 
value 9000 (anything >4600 
will work, but 9000 is a 
quick and easy number to 
enter > . 

II. If it is desired to raise 

the probe, enter the extreme 
value 90 (anything less than 
940 will work ) . 

e. -The program will next ask which bit 

IS selected to be "high". The 
following guidelines apply: 

i. If it IS desired to lower the 
probe, enter 7. 

ii. If it IS desired to raise the 
probe, enter 3. 

f. Switch ON probe activation power and 
when the probe reaches a desired' 
degree of deflection, switch the 
activation power OFF and record the 
mV value printed on the screen. 

g. Repeat steps (d) through (f) until 
sufficient data is collected. Load 
this data into the program "POLYFII" 
and request a second order fit 

(let X=angle and Y=mV). 

Enter the coefficients derived by 
this program into the appropriate 
location in the program "PROBE_SUBS" , 



NOTE: The relays will be activated on 



88 



610 

620 

630 

6^0 

650 

660 

670 

630 

690 

700 

710 

720 

730 

740 

750 

760 

770 

780 

790 

800 

810 

820 

830 

840 

850 

860 

870 

880 

890 

900 

910 

920 

930 

940 

950 

960 

970 

980 

990 

1000 

1010 

1020 

1030 



? digital 1o(j! When the machine boots 

? up (hp-9826), all relay lines are 

? high (+5V). The instructions belou 

1 uill drop the voltage to zero. 

» 

INPUT ’*What is the desired position? ( mV ) ’* . Vo 1 tage 
PRINT "Desired position (mV) *";Voltage 

t 

IRead the actual motor position. 

? IF the desired position is BELOW the 
? actual position, then tell the motor 
! to move UP. 

? If it is ABOVE, then tell it to go DOWN. 



INPUT "WHICH BIT DO YOU WANT HIGH7",Lbit 



OUTPUT 723;"0P,7,0T" ?CLEAR ALL D/0 
OUTPUT 723;"0P,7,";2^Lbit;"T" 

t 

OUTPUT 723;"CCJT" ?CLEAR A/D CARD 

OUTPUT 723;"SF,3,3,3. 1 .25, 12T" !F0RMAT A/D 
OUTPUT 723;"0B,1 ,10,1T" ?CL0SE THE RELAY 

!THAT CONNECTS THE 
!A/D CARD TO THE 
♦POTENTIOMETER 

OUTPUT 723:"IP,3T" ?START A/D CONVERSION 
ENTER 72301 ;A lENTER A/D VALUE INTO A. 

DISP A 

BEEP ABS(A) . .01 

f 

? The follouing IF stops the motor uihen 
I it reaches it's limits. 

t 

IF (Lbit-7 AND A>4600) OR (Lbit-8 AND A<940) THEN 
OUTPUT 723;"OP,7,0T" 

BEEP 

PRINT "ALL DONE??!" 

END IF 
GOTO 890 
END 
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7 . MOTOR CAL 



10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
no 
120 
130 
140 
150 
160 
170 
130 
190 
200 
210 
220 
230 
240 
250 
260 
270 
-) ) 
280 
290 
300 
310 



fiOTOR^CAL 

PURPOSE: Calibrate the potentiometers 

used uith the motors on the milling 
mach 1 ne . 



OPTION BASE 1 
DIN Coef<12) 

f 

OUTPUT 723;’‘SF.3,3,3n .25. 12T’* 
LOADSUB ALL FROM ••MTR_SUBS” 
LOADSUB ALL FROM ’’PROBE^SUBS*’ 

F 1 1 eS = "mo tor_coef '* 

I 

CALL Cal ibrate<Fi le$> 

CALL Retr leve^coef <Coef (*) ,Fi leS) 
BEEP 2400 J 



PRINT ’*What (x.y.z) position do you wish to’* 

PRINT ” move to? (inches relative to noz-” 

PRINT •* zle)” 

INPUT X,Y,2 

f 

CALL Move_ldv_to (X , Y ♦ Z , Leng th_arm , Length_probe , Angle_arm , Probe_ang le , Coe 

f 

BEEP 

GOTO 280 
END 
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8 



LOAD XYZ 



10 


1 


lGAD_XYZ 


20 


f 




20 


1 


by Bill Calbreth 


40 


1 




50 


f 


30 April 1984 


60 


1 




70 


f 




30 


1 


PURPOSE: ’This prograjD uiill ailoui the 


90 


i 


user to enter desired <x.y,z) posi- 


100 


1 


tions of the milling machine relative 


110 


» 


to the nozzle. The values will be 


120 


j 


stored on disk to be utilized later 


130 


\ 


by MAIN_T. 


140 


\ 




150 


DIM X<500) . Y(500) ,Z(500) 


160 


I 




170 


» 


1. Input the file najr>e . 


130 


1 




190 




GQSUB Clear_screen 


200 


j 




210 




PRINTER IS 1 


220 




PRINT "Input the file najDe." 


230 




PRINT "{ I suggest 'RUNXX' uihere 'XX'"^ 


240 




PRINT " is the run number. 1" 


250 




PRINT 


260 




INPUT F i lenameS 


270 




f 


280 


1 




290 


1 


2. Begin inputting data. 


300 


I 




310 


1 




320 




GOSUB Clear screen 


330 




1 


340 




PRINT "Do you Wish to append a previous ( 


350 




INPUT Answers 


360 




1 


370 




IF AnswerS-"YES" THEN 


380 




INPUT "Previous file name'?", Old fileS 


390 




ASSIGN 5kFilel TO Old fileS 


400 




1 


410 




1 = 0 


420 


Loop! : ! 


430 




I = T>1 


440 




ENTER ^Filel ; X ( I ) , Y ( I ) . Z ( I ) 


450 




OUTPUT 701 ;"I ,X,Y,Z = ";I,X<I) .Y(I) .Z(I 


460 




IF X<I)O-100 THEN GOTO Loopl 


470 




ASSIGN ivFile] TO * 


480 




PURGE 01d_fileS 


490 




Coun t = I - 1 


500 




ELSE 


510 




Coun t = 0 


520 




END IF 


530 




1 


540 




GOSUB Clear screen 


550 




BEEP 


560 




» 


570 




PRINT "1. Input the desired position in 


580 




PRINT " inches as X . Y, and Z." 


590 




PRINT 


600 




PRItir "2. Terminate input by entering 
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810 
620 
630 
6^0 
650 
660 

670 Begin: ? 
680 
690 
700 



PRINT 

PRINT 

PRINT 

PRINT 

PRINT 

f 



Terminate input by entering '-100'” 
for X , Y . and Z . “ 

If you ijiish to enter an orientation angle,'* 
enter ' -999 , or len tat ion , 0 ' . ” 



Count^Coun t+1 
PRINT ‘’Item <f’*;Count 

INPUT '*(X,Y,Z) in inches?” , X (Count ), Y (Count ), Z<Count ) 
OUTPUT 701 ; ”I.X,Y.Z-*’;Count,X( Count ) .Y(Count) ,Z(Count) 
IF X ( Count )<>- 1 00 THEN GOTO Begin 



End_count*Count- 1 

All data points have been entered. 



a. 

b. 



GQSUB Clear^screen 

f 

PRINT ’’SOFTKEY LABELS; 



710 
720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 
960 
970 
980 
990 
1000 
1010 
1020 

1030 Clear_screen : 

1040 
1050 
1060 

1070 OUTPUT 2 USING 

1080 GCLEAR 

1090 RETURN 
1100 
1110 
1120 
1130 
1140 
1150 

1 1 60 Hard_copy ; 

1170 ! 

1180 
1190 
1200 



Set up new softkeys. 
Explain sof tk eys . 



WRITE data to disk.” 

EDIT out bad data.” 

HARD copy the data on printer,” 
LIST data on the CRT.” 

STOP terminates the program.” 



PRINT 
PRINT 
PRINT 
PRINT 
PRINT 

t 

ON KEY 0 LABEL '’WRITE” GOSUB Write.data 
ON KEY 2 LABEL ’’EDIT” GOSUB Edit^data 

ON KEY 4 LABEL ’’HARD” GOSUB Hard_copy 

ON KEY 6 LABEL ’’LIST” GOSUB List_data 

ON KEY 8 LABEL ’’STOP” GOTO Terminate 

GOTO 900 



Clear the CRT display. 

-s^,B”;255.75 



! Print out all data to the printer. 
PRINTER IS 701 
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Y C coan c ) 



Z < count ) 



1210 ! 

1220 PRINT '* Count X< inches) 

1230 PRINT •• 

1240 PRINT 

1250 FOR 1=1 TO End_count 

1260 PRINT I.X(I) ,Y(I) ,Z(I) 

1270 NEXT I 

1280 ! 

1290 PRINTER IS 1 

1300 RETURN 
1310 ! 

1320 ! 

1330 ! 

1340 ! 

1350 ! 

1360 ! 

1370 List_data: ! 

1380 ! 

1390 ! List data to the CRT. 

1400 ! 

1410 GOSUB Clear_screen 

1420 » 

1430 PRINT “There are" ; End_coun t ; data points.'* 
1440 PRINT 

1450 INPUT “Which point do I start u i th?“ , S tar t 

1460 INPUT “which point do I end u i th?“ , Ena_data 

1470 ! 

1480 FOR I-Start TO End_data 

1490 PRINT “I.X.Y,Z = “;I,X(I) ,Y(I) ,Z(I) 

1500 NEXT I 

1510 RETURN 

1520 ! 

1530 ! 

1540 ! 

1550 f 
1560 ! 

1570 ! 

1580 Edit_data; ! 

1590 ! 

1600 ! Edit data. 

1610 ? 



1620 


GOSUB 


Clear_screen 






1630 


1 








1640 


INPUT 


**Which data point 


do you 


wish to ai ter'?*' , I 


1650 


» 






1660 


PRINT 








1670 


PRINT 


“For point ^“;I;“ 


. (X.Y. 


Z) where:'* 


1680 


PRINT 


X<I) ;",“;Y'(I) ;**,“ 


;Z<I) 




1690 


» 








1700 


INPUT 


“Type in the new 


values : 


**,X<I) ,Y<I) ,Z(I) 


1710 


RETURN 






1720 


» 








1730 


f 








1740 


t 








1750 


» 








1760 


1 








1770 


\ 








1780 


Write data: 


♦ 






1790 


» 








1800 


! Wr 


1 te data out to a 


file on 


d 1 sW . 
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1810 

1320 

1830 

18^0 

1850 

1860 

1870 

1880 

1890 

1900 

1910 

1920 

1930 

19^0 

1950 

1960 

1970 

1980 

1990 

2000 

2010 

2020 

2030 

20^0 

2050 

2060 

2070 

2080 

2090 

2100 

2110 



X < End_coun t + 1 )--100 
Y ( End_coun t + 1 ) * - 1 00 
Z ( End_coan t + 1 ) 1 00 

I 

Maxidata* 3* < End_coiin t”*- 1 > 

» 

F 1 le_s i ze“INT < Max_data*8/256 ) ■►1 
\ 

CREATE BOAT F i lenameS ,Fi le_si ze 
ASSIGN «File TO FilenaineS 

FOR 1*1 TO End_coant'*'l 

OUTPUT @FLle;X(I) . Y(I) ,Z(I) 

NEXT I 

» 

GOSUB Clear_screen 
BEEP 2^00, .3 

PRINT "File " ;F i lenameS ; " has been stored 
ASSIGN aFile TO * 

RETURN 



Terminate: ? 

GOSUB Clear_screen 

PRINT "NORMAL TERMINATION OF PROGRAM? 
END 
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SEND DATA 



10 
20 
30 
40 
50 
60 
70 
80 
90 
100 
1 10 
120 
130 
140 
150 
160 
170 
180 
190 
200 
210 
220 
230 
240 
250 
260 
270 
280 
290 
300 
310 
320 
330 
340 
350 
360 
370 
380 
390 
400 
410 
420 
430 
440 
450 
460 
470 
480 
490 
500 
510 
520 
530 
540 
550 
560 
570 ^ 
580 
590 
600 



SEND.DATA 

To VAX.IBM,TRS-80. 

HP-9826 TERMINAL PROGRAM 
[REQUIRES BINARY ENHANCEMENT PROGRAM 
”BE3'’?] 

JUNE 30. 1982 
updated 1/5/83 
updated 1/16/84 

BILL CULBRETH 



Sc*9 ! RS-232 IS SELECT CODE 9. 

PRINTER IS 1 ! PRINTER IS CRT. 

Pr-1 ! DEFAULT PRINTER IS CRT 

Pr inter_choice*701 ! MY PRINTER IS 701. 

Bits-7 ! BITS PER CHARACTER 

Duplex-0 ! FULL DUPLEX 

Baud-300 ! BAUD RATE 

Computer-1 ! ASSUME IBM COMPUTER 

I 

OUTPUT Pr;*'{300 BAUD. IBM assumed.” 

OUTPUT Pr:” Load the binary program BEB first” 
OUTPUT Pr;” unless you have BASIC 2.0” 

OUTPUT Pr;” SET MODEM ON <FULL DUPLEX> }” 

OUTPUT Pr;” ” 

f 

DIM Name$C200] .Hp_file$C 301 .Aa< 1500) ,NumbS(301 
INTEGER Isend 

f 

CONTROL Sp.3;3aud 

CONTROL Sc.4;Bl ts-5^^4 ! BITS/CHAR A ^^STOP BITS. 

» 

f 

to_d isk -0 

Datadump-0 

I_data= 1 

1 = 1 

J=1 

K=1 

L-1 

ON ERROR GOTO Errors 
ON KEY 0 LABEL ’’Line Mode 

5 LABEL ”Termina 

6 LABEL 

7 LABEL 



ON KEY 
ON KEY 
ON KEY 
ON KEY 8 LABEL 



GOTO Line_mode 
GOTO Terminal 
'To Crt” GOTO Pr_crt 
'To Prt” GOTO Pr_prt 
'DATA” GOTO Data_dump 



Line_mode: ! 

OUTPUT Pr;”{LINE RECEPTION MODE}” 
Begin; STATUS Sc.l0;Y ! CHECK FOR FULL BUFFER 

IF BIT(Y,0)-0 THEN GOTO Begin 

I 

! RECEIVE ROUTINE. 
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610 ! 

620 Receive: 



STATUS Sc,6;A 



630 
640 
650 
660 
670 
680 
690 
700 
710 ^ 

720 
730 
740 
750 
760 
770 
780 
790 
800 
810 
820 
830 
840 
850 
860 
870 
880 
890 
900 
910 
920 
930 
940 
950 

960 Data_dujnp : 

970 
980 
990 
1000 
1010 
1020 
1030 
1040 
1050 
1060 
1070 
1080 
1090 
1100 

1110 GOTO Begin 
1120 I 
1130 ! 

1140 ! ERROR HANDLING SUBROUTINE 
1150 ! 



B= A 

OUTPUT Pr USING , A” :CHR$(B) 

IF B=*63 AND Datadojnp*1 THEN GOTO Data_dojnp 
IF B-13 AND Computer>3 THEN OUTPUT Pr:CHR$<13) 
GOTO Begin 

TRANSMIT ROUTINE. 



IF Duplex-0 THEN 

IF NUM<Key$)<>255 THEN OUTPUT Pr USING ’V^A” 
IF NUM<KeyS)>255 THEN OUTPUT Pr;” ” 

END IF 

IF Computer-1 AND NUM(Key$)-8 THEN Key$-CHR$ ( 64 ) 

the previous line gives an ^ 
for a baclfspace for the IBM. 

IF Computer-5 AND NUM<Key$)-8 THEN Key$-CHRS( 1 27) 



THE UAX/VMS REQUIRES A DELETE 
SYMBOL FOR A BACKSPACE. 

IF NUM<Key$)-255 THEN Key$-CHR$( 1 3 ) 
OUTPUT Sc USING ”#,A”;KeyS 
GOTO Begin 



DATA FILE OUT TO THE HOST COMPUTER. 



IF I_data-1 THEN GOSUB Open_fiie 

f 

IF Datadump-0 THEN GOTO Begin 
IF Computer-1 THEN WAIT .3 
! wait for the slow IBM. 

BEEP 1000+RND*1500, .05 
OUTPUT Pr;”A(”;I_data;”)-”; 
OUTPUT Pr ;Aa(I_data) 

GOSUB Send_number 
IF Aa(I_data)--200 THEN 
I_data- 1 
Datadump-0 

END IF 

I_data-I_data+1 



1160 Errors : 


OFF ERROR 


1170 


Close f i 1 e-”200 


1180 


! FIRST. END OF FILE 


1190 


IF ERRN-59 THEN 


1200 


Aa< D--200 



;Key$ 
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1210 
1220 
1230 
12^0 
1250 
1260 
1270 
1280 
1290 
1300 
131 0 
1320 
1330 
13^0 
1350 
1360 
1370 
1380 
1390 
MOO 
MIO 
M20 
M30 
1440 
1450 
1460 
1470 
1480 
1490 
1500 
1510 
1520 
1530 
1540 
1550 
1560 
1570 
1580 
1590 
1600 
1610 
1620 
1630 
1640 
1650 
1660 
1670 
1680 
1690 
1700 
1710 
1720 
1730 
1740 
1750 
1760 
1770 
1780 
1790 
1800 



GOTO 2000 ! RETURN AFTER ERROR. 
END IF 



IF ERRN059 THEN OUTPUT Pr;’*<error ^’‘lERRN;’* generated . >” 

IF ERRN = 54 THEN OUTPUT Pr;‘*{FILE <“ ; Hp_f i ie$ ; “> ALREADY THERE)" 
IF ERRN-54 THEN GOTO Created 

IF ERRN=56 THEN OUTPUT Pr;"<FILE <" ;Hp_f i le$ ; " IS NOT ON DISK,)" 
ASSIGN ^File TO * 

GOTO Line_mode 



i OUTPUT TO CRT. 



Pr_crt: Pr-1 

GOTO Line_mode 

f 

I 

! OUTPUT TO PRINTER. 

I 

Pr_prt: Pr-Pr in ter_cho ice 

GOTO Line_mode 

t 

f 

! CHANGE THE TERMINAL CHARACTERISTICS. 

f 



Terminal : ! 



OUTPUT 


Pr ;" 


1 . 


Baud Rate -";Baud 


OUTPUT 


Pr;" 


2. 


Bits/Char =":Bits 


OUTPUT 


Pr;" 


3. 


Duplex -"‘.Duplex 


OUTPUT 


Pr ;" 




[ 1 -f ul 1 ,0-half 1" 


OUTPUT 


Pr :" 


4. 


Computer -";Computer 


OUTPUT 


Pr;" 




[IBM-1 . UAX/UNIX-2. " 



OUTPUT Pr;" TRS-80-3. Cyber-4, s/ax/vms-53" 

OUTPUT Pr;" " 

INPUT "Change which one7", Which 
IF Which“l THEN INPUT "To9”.Baad 
IF Which=2 THEN INPUT "To?". Bits 
IF Which-3 THEN INPUT "To?" . Dup lex 
IF Which-4 THEN INPUT "To?" . Computer 
IF Computer*! THEN Duplex*0 
IF Computer*3 THEN Duplex*0 
IF Computer-3 THEN Bits-8 
IF Computer-5 THEN Duplex-1 
GOTO Line mode 



t 

j 

Open_f i le : ! 

f 

I 

f 

Datadump- 1 
f 



Open a file to read data from 
disk . 



INPUT "Is this LDU data? { 1 - YES) " . Ldv$ 

IF Ldv$-"1 " THEN 

INPUT "Experiment . Exper imentS 
ELSE 

OUTPUT Pr;"Data file out of HP to host." 
INPUT "File name'^" , Hp_f i 1 eS 
END IF 
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IBIO ! 

1920 IF Ldv$»’M“ THEN 
1830 Hp_f 1 le$*Exper imentS6’'_RESULT” 

18^0 END IF 

1850 ! 

1860 ! Read the file off of disk. 

1870 ! 

1880 ASSIGN ^File TO Hp_file$ 

1890 I»1 

1900 Check-0 

1910 BEEP 

1920 BEEP 

1930 OUTPUT Pr ;”(Working on file <*’ ;Hp_f i le$ ; 

19^0 ! 

1950** 

1960 ENTER @File:Aa<I) 

1970 Check-Aa<I) 

1980 I-I^-l 

1990** 

2000 ? 

2010 ASSIGN i^File TO * 

2020 Datadujnp-1 

2030 RETURN 
2040 ! 

2050 ! 

2060 ! 



2070 Send namber: ? 



2080 

2090 

2100 

2110 

2120 

2130 

2140 

2150 

2160 

2170 

2180 

2190 

2200 

2210 

2220 

2230 

2240 

2250 

2260 

2270 

2280 

2290 



RETURN 

I 

I 

t 

END 



t SEND A NUMBER ONE CHARACTER AT 
! A TIME TO THE HOST COMPUTER. 

f 

NiuT>b$-UAL$( Aa< I^data) ) 

Lengt h=LEN< NujnbS) 

I 

IF <<Ldv$-"r’> AND <I^data>13)> THEN 
Posit-POS<Niunb$.”.”) 

IF (PositOO) THEN Length-Posi t+2 
END IF 

I 

FOR 1-1 TO Length 

Niuneric-NUM(Najnb$I I .11) 

OUTPUT Sc USING , A” ; NiunbSI I , I ] 
NEXT I 

I 

OUTPUT Sc USING , A” ;CHR$< 1 3) 
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APPENDIX C 



MAINFRAME PROGRAMS 



1 . TCAL 



tcal 

purpose: this POCGNAM PLaCcS ThERMOCCUPLE calibration data 

RtCEl VFO FROM the HP-<3826 MICROCOMPUTER INTO A 
MORE workable format for USE IN THE PROGRAM TF I T • 



OI mens I ON TPROBEI 200l.rAMe(200l .r NO 22(2 00 I 
I = I 

10 RE AD( 0 7 . $ ) r AMR( M 

IF( TamB ( I I .EQ.- 200 . 0 I GO To 20 
REAOf 07 I TPRObF I I I 
RE aO( 07 ,$ ) TNG22 ( I I 
NI fEMS= I 
I = I ♦ I 
GO TO 10 
20 CONTINUE 

A w I T E ( N . A O I 
DO 30 I = 1 .N I TEMS 

^»RITECu.SO» I .TAmhI I I.TPOOBEI I I ,TNn2Z(I I 
30 CON r I NUE 

AO FQPma r I bX . • AM0 I tN T I C I • . 2X . ‘PROBE ( C I * . 2 X , • N02 ZL £ (CMI 
SO forma r ( I X . I 3 . ox , 2 I P6 . J . oX M 
STOP 
ENO 
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2. TFIT 



rFi r 

PUWPQSf: THIS PQOGO^M PERFORMS a FIRST UROER CURVE Fir BY ThF 

LEAST squares method fqq r»c cclumns of data. the 

FIRST column lists VALUES OF X anO TmE SECOND LISTS 
values of F(XI or Y. for THER»-0C0UPLE CALIRRATIUN, 
let X = THt AM0itNT rFMPeoArcRt* and let ^ - either 
THE PROBE OR NQ^^lE TEMPERATURE aS DESIRED. 

dimension x| I SO I .Yl 1501 .XSO( I50I.XYI 150 I.YESTI 1501 

SUMX =0 . 0 

SUMxSOsO . O 

SUM XY = O . 0 

Sumy = 0.0 

SUMnuMs 0 • O 

SUMDE N= 0 • U 

RE ADI 7 . I N I TEMS 

DO 10 I = I ,N I TEMS 

WE AD ( 7 , c: I X I I I • V I I I 
XSQI I » = X ( I I 
X V ( I I s X ( II c: Y U I 
Su MX =SOM X ♦ X I I I 
SUMXSO aSUMxSQ ► xsc ( I » 

SUmxy = SUmxy*-xY( I I 
SUM Y = Sum Y ♦ Y I I I 
10 CONTINUE 

XI TEmssFLUATINI TFMS I 

Ar( SUMYiSSUMxSC-SUMX^SUMXY |/IX lTeMS-*?SUVXSQ-SUMx^*2| 

• 3 ~ { XI T'^MS^SIJMxY-SUMX^SUmy |/( X (TEmS=>SUMXS0-Sumx<:$2 » 
yijarzsumy/x I rt^s 
DU 4.>0 I s I • N I TE«S 

YE5 r I I I =rf* X ( I I ♦ A 
#RlTF(6.^:»VESrlM. v(ll 

SUMNUMsSUMNUM * (VEST (I l-yHAW)^t:2 
SUMOE NZSUMUEN«-IY( II-YBA«|^::?2 
20 CON r I NUE 

RSO = SUMNUm/ SUMOE 
«R I TE I O • ••:: I A ,0 ,RSQ 
S TOP 
END 
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3 . JETCURV 



J6TCURV 



PURPOSE: THIS program fits jet TRAJECTCRY data to the equation 

Z=A$Y/(U**Y| 8Y the least SOUAWeS mETHOO and computes a 
COR wOLATION coefficient DATA POINTS APE READ 

IN free format from unit "7” C I SK , 



dimension YI50I*ZI50).yy( 50),ZZ(50I •YSQ(50).YZ(50)»ZESri50) 

SUMYsO . O 

SumySO=0 .0 

SUMYZ =0 .0 

SUMZ =0 , 0 

SUMnU M = 0 . 0 

SUM0EN = 0 . 0 

READ! 7,4:)Nl TEMS.FACrCR 
DO 10 l=l.NIT£MS 

READ! 7,:S»Y(1 ) , Z{ I \ 

Y I n =Y ( ( )/F actor 

Y Y I I 1 = 1, 0/Y I ( ) 

ZZ( I »=l.U/IZ(I I/FACTORI 
YSO ( I ) -YY ( I 2 

YZ I I » = YY (I ) $ZZ I I ) 

SUMY = SUMY* Y Y ( M 
SUMysQ = SUMYSO •> YSU I I ) 

SUMYZ =SUMYZ ♦■YZ ( I \ 

SUMZ =SUMZ ZZ I I i 
10 CONTINUE 

KlTcMS=FLGATlN( TENS) 

AA=(SUMZ^:SUMYbQ-SUMY^:SUMYZ)/( XlTEMSC:SLMYSO-SUMYi:^:2 ) 
t3b = I X I TE M S 4* SUMYZ- SUM Y:?S0MZ |/( XireMS*SLMYSO-SUMY^:02 ) 

A- I . 0/ A A 

e = 

ZHAW=bUMZ/X ITEMS 
00 20 I - I . N [ TEMS 

ZES Till =EOi:YY I n ♦ A A 
^RITtlo.JOJZZI II. ZEST! I) 

SUMNUM = SUMNUM «‘IZESr I l)-ZbAR|*0 2 
SUMOLN-SUMDEn* (ZZI I )-ZBAWI^:52 
20 CONTINUE 

RSO=SumnUM/SUMOEN 
<W I TE I 6 . AO » A , U . RSC 

JO FURMAflix.'Z (actual! = ' . F I 2 . 8 . J X . • Z (EST) =*.Fl2.8! 

AO FORM A r I I X , • A= • .F 12 .8 . 2X . • . F I 2 . d . 2 X , « R 2- * .F I 2 . d ) 

s f np 

END 
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GRAB 



c 

C grab 

c 

c purpose: data transfer from the HP-Reze'^ i crocqmputor 

C TO THt IBM, 

c 

C BY BILL CULdPtTM 

C FOR ^EZaIO. fall GUaPTER. 1982 

c 

C F ilEOEF 05 terminal 

C FILEOEF 06 terminal 

C FILEOEF 07 DISK mvqaTa DATA ( PERM | 

C 

C global TXTLI0 F0RTM002 MOC2EEH 

C 

C TYPE IN THE ABOVE LINES TO MAKE THIS 

C FORTRAN program Run. 

c 

c 

dimension oataijoooi 
I = I 

V»R I TF ( 6. 80 I 

80 FORM A I ( 2X , • Q6C I N INPUriNC Data From the hp-9026*I 

C 

C 

10 CONTINUE 

RE AUI 5 • C: I OAT A M I 
I = I ♦ I 

IF ( DA TA ( I- M .nE .-200 » GOTO lO 
C 

N I TEMS - I - I 

6 FORM A r ( 2X . I 5 . * DATA POINTS xERE eNTEREO.*! 

mR I IE I 6. 6 » N I TEMS 
C 

C NOW THAT ALL DATA HAS BEEN ENTERED. mRITF IT OUT ON 

C DISK. 

C 

S FO3MAr(2x,»0ATA(*,l5.M = ••IF15.5I 

I = I 

20 -R I TE I 7 I OAT A ( I I 

I = I * I 

I F I OA r A ( I - M .nF .-200 I GOTC 20 
C 

C ALL DATA HAS BEEN /xRITTFN ONTO DISK. 

C 

STOP 

END 
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5 . TDATA 



rc A r * 

PUOk^GSE: Th ( s PC?CGSA>- PLACES ThE =3<jC>ant JET OATi 5 

F^^n*- r^E Hp-qb2e « i c»ccc'-pu te i stc a *-upe 

FCR«AT. IT alSC CCnvERTS U^lT5 CP LEnGTh 
TO •'- E r » I c . 

01 “ENSIGN XI200I.Y[200I .ZI200I.rP»ObEI200I .rA*-0i2OOI 
ISTOEVI200I .XM«I20C» .rM«i200» .Z“-( 2QQ» 

( = I 

10 REAOI 07 ,ci X II I 

IFI Al M .EG. -200.0 » GC TO 20 
I I 1=25 I I I 
RE AU( o 7 ,c) y ( ( » 

Y*-« ( I I = 25 . A^ V ( I I 

REAUI07.:5I ZIM 
Z““ I I 1=15. AcZ I ( I 
RE AOI 07 , I TDROeE I I I 
R£AC( C 7 I T A-H I I I 
READI oz .= I TNQZZI I I 
REACI 0 7 I 5TCEV I I I 
N ( TE“S= ( 

I = I ► I 

GO TO 10 

20 continue 

• N I TE I 5 . aq I 
00 JO I = I . N I TE“S 

-RITEIo.SCI f .<«“! n .YMwf I I.Z-wlI I .TORGbEl I ). TAugi 
t STOcV ( I ) 

JO CGNTlNUt 

4 0 PGR -'Aria*,*x IWM|*,4X,*Y (•‘•“l•,4X,•z l““I*.2*.* 

I'AWblENT ( Cl • . I X , » NCZZLfz (CI'.lx.'STO GcV (Cl*» 

SO c^GRVAri IX . IJ.Jx.2lP7.J.JX».»- 7.j,ax,^(F^S.J,AX | | 

S TCP 
ENO 



ECEl VEO 
CPOFPL T 
FRG“ 51 

. TSOZZ 12 00 I . 



I I . TnCZZ ( M . 

IC I • . I X . 
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6 . C 0 NT 0 UR 4 



C 

C 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 

c 



c 

c 

c 

c 

c 



c 

c 



c 



c 

c 

c 

c 

C^: 

c 

c 



c 

c 

c 

c* 

c 

c 



CON 



PURPOSE: THIS PRCGPA^* is OESIGNEO TC C I SPL A Y BUOYANT jFT 

TE '-OfcR A TUPE DATA USING A CONTOUR PL Q T T I NO PACKAGE 
avaIlaelE ^ITH OISSPLA. THE FOLLOWING Raw OaTa IS 

read FRCM a 0 I S< file: PROEE positions in xyz 
COORUINATES relative to the nozzle tip, AMhIENT 
TEMPERATURE. PROBE TEMPERATURE AnO NOZZLE IFmPERATURE. 
THE PRQPE ano nozzle T Em PE P A fUwE S aRE TRanSFORmpu 
BY calibration coefficients anu normalizeo ww t ambient 
temperature. the xyz CQURCINATES are converteo to 
xSw COQROInateS RFlaTIVE TC the intersection of the 
data plane with the centerline trajectory uf thf jet 
(The S-axIS is tangent to the TRajECTOWYI. this 
program also CCmPUTES the raie of hfat transfer from 
THE JET TO THE ambient. 



note: the following values must be inserted as the 

FIRST line of CATA in free FCR-AT: ThF total number 
OF DATA POINTS INITEMSI. THE ACUTE AnGlE RETwEEN 

The OAfA Plane and horizontal (TmETai ano the vertical 
OISTANCE BET.IEEN ThE Z-axIS ANO THE INTERSECTION 
POINT CISCUSSEC ABOVE (ZA|, THE CENTERLINE VELOCITY 
IN M/s (VELl ano jet width in mm (WIOTHI, 



dimension XIlOOl.YllOOl.TIIOOl.TPIlOOl.TAllOOl.TNI IOOI.TM aTI lO.lOl 
common wORKllbOOOl 



* c: READ DATA FROM UNIT '* 7 •• DISK 0^:0: 



REAOC 7.^1 NlTEMS, THETA.ZA.VEL.wICTH 



00 20 I=I 


.N I TEMS 






REAOI 7 


. * 1 A ,X 1 


( 1 • A • Y I 11 


. TP C ( 1 , T A ( I 1 . TN 1 I 1 • A 


Y I I 


1 - ( Y ( I 1 


-ZA I/SIM ThETAI 






calibration cnEFFICIENTS 


TP( 


I 1 = . RBAEORe 75<:rp I I 


1 1 . 74 07qt)C5 


TNI 


I 1 = 332234 2e3=^TN I I 


1 1 . d4 75 40 66 


T 1 I 


1 = ( TP I I 


1 - ta I I 1 


1 /( TNI I 1 - TaI I 1 1 


20 continue 








30 format I 2X 


t *I.X,Y,T 


- • • IS . 3F 1 


5. 51 



DIMENSION OF THAT ^ $ $ 



I XO I M = 10 

I YO I M - 10 



•4; :> 4-' FIND The MAXIMvjM ano minimum TEMpFRATURES ^ ^ ^ 



T M I N = I 00 . 0 
TM AX = 0.0 

XMAXsQ. 0 



YMAX=0-0 
XM I Nr I 00 . 0 
YM I Nr I 00 . 0 

00 40 I r I , N I TEMS 

IF(T(II .CT. TMAXI TMAX = Till 

IFiTin .LT. rviNl Tmin s TIM 
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40 

c 

* 

C 

C 

C 



C 

C- 4 
C 

C 

c* * 
c 

c 

c 

c* $ 

c 

c 

c 



c 



c 

c 

c 



c 

C;:: * 

c 



c 

c* * 
c 



(F(X(I »,GT.XyAX|XMAX = X(l ) 

IF(Y( I I ,GT.VMAX|YMAX = Y( I I 
lFiX(I».LT.XMlN)XM|N:=X(() 

|F(Y( ( ) .LT,YM(MYM(N = Y| I \ 

CON r I NUc 

^ $ OETEH*^lNk£ the maximum VALUES OF x anu Y •#ITm(N A 

OUAORANf, ASSUMING The STREamm(SE AXIS IS CENTEf^EO 
IN A SQUARE matrix Plane 

0x = ( (XMAX-XMInI/FlOATI IXOIMI I/IQOO. 

0Y = I (ymax-ym(n»/FLOAT| IyDIM) )/ioOQ, 

QUA0X=(FLUATI lXUlM|/2.0)^iCX 
QUAO Y = I FLOA r I I YO I M I /2 . 0 » ^CY 

* find the incremental area for heat transfer calculations 

ARE A=OXCO Y 

* establish the increment SIZE FQR CONTOUR PLOTS * C: ^ $ 
TINCR=( Tmax-TmINI/5.25 



* ^ generate the grid $ * $ 



call COMPRS 
Call tek6 i e 
Call page ( d. . s . 1 

CALL BLU«uP (0,751 
CALL PHY SQR 11.25.1.1 
CALL AREA20 I 5 . 5 .5. 5 » 

CALL HE I GH T I 0 . 100 1 
CALL CARTUG 
CALL INTAXS 

CALL xnameI'X (mm|s»,I001 
CALL YNAMEI'W I MM I $ • , I 00 I 

C all X T ICKS ( 2 » 

C all y T I o< S I 2 I 
call X a X anG I «50 . 0 I 



CALL SaISSM 

CALL ME AO I N I • temper A TLRE CONTOURS IN A EUUYANT jETI*. 100.2. 
Call hFAOIN|»«ITh a CHOSSFlQwinG amb I EN T t • • I 00 . 2 . . 4 » 

CALL hEaOIn(»(25X FlCa R a tE ) -S • . I 0 C . I . . 4 » 

CALL mE AO ( N { • I PL ANE A ) t | 0 0 • 1 • • 4 | 

CALL GRAFI-3Q, 0.05.0. JO. 0. -JO. 0.05. 0.20. 01 

* GENERATE The InTERPOLATEO TE mo ER a T LP j:: maTR(X "TmaT” ^ * 



CALL bCOMUN ( I oOOn ) 

call zb ASE ( TM I N I 

call PGNm a t I I xo I m , ( yo I m ) 
call GE Tm a r ( X . Y , T . N I TE ms . 0 ) 
call E noma T ( TM A T , 0 I 

* * SOLVE the rate of heat TRANSFER * $ i ^ 



I 

1 



QSUM=0 . 0 

UQ 70 I=I.IXUIM 

00 bO J = I . I YO I M 

Rn=ll.00004jj . J.SlRlOE-SbTMATII.Jl 

-5.Qb^Zbo7E-t^TMAT I I ,J )bC2 
J. 4 I 2 ')o 29E-P<-. T M A r ( I , j )c-^ J 
-N.740b7LF-l 14:T4 aT ( I .J)*^4)$|0 0O. 
XX^QUAUX-FLOATI J-I |-::0X 
Y Y=OUAr.Y -F lO A r ( (-1 J?:LY 
R AO I US =SLR T I X XC -? 2 ♦ YY<:b2» 
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non 



V=vEL^&XPl-PAOlUS**2/v«lOTHee2 1 
(J = w0^ AWE A^VC.4<*5^(rMAT|i .J»«.273.16» 
OSUMrOSUM ♦Q 
60 CON r I NUt 

70 CONTINUE 

i-P I TE ( <5 • tiO I OSUH 
80 FOPMAT ( I X . '0= '.F20.5J 

^ ^ ^ PLOT The CQNTOUPS ^ ^ 

CALL C0NMAK(TMAT, IX01M,IYDIM,TINCP» 

CALL CONLiN(0.*SOLlCS* .'LARELS* .1 .101 
Call COnanG(po.oI 
CALL PASPLN(0.2b» 

Call CON tup ( i .* labels •. ‘Dpaw • I 
CALL height (0.05J 

CALL WLMEbS(*EXP IC t • . I 00 .-59. -59 I 

CALL DOT 

call CP I 0( 1.11 

CALL LNOPL(0» 

CALL UONEPL 

STOP 

ENO 
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10 

15 

20 

25 

30 

35 

40 

45 

50 

55 

60 

65 

70 

75 



APPENDIX D 



TABULATED DATA 
TABLE 1 

ROTOMETER CALIBRATION 



ml/s 


Std Dev (ml/s) 


5.94 


0.32 


8.02 


0.56 


10.05 


0.41 


11.85 


0.24 


14.03 


0.33 


15.98 


0.58 


17.41 


0.31 


19.09 


0.34 


20.75 


0.23 


22.38 


0.28 


24.27 


0.22 


26.03 


0.20 


27.64 


0.37 


29.53 


0.24 
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TABLE 2 



TEST RESULTS 



Crossflow velocity: 




.13 m/s 


Nozzle 


flow rate: 




11.85 ml/s (25 


Nozzle 


inside diameter: 




7.144 mm 


Nozzle 


discharge velocity 


(mean) : 


29.558 mm/ s 


Nozzle 


temperature (mean) : 




41.8 °C 


Ambient 


temperature (mean) 


: 


24.9 °C 


Fro ude 


number : 




14.8 


Michaelis-Menter equation: 








Coefficient A: 




2.64284325 




Coefficient B: 




1.03698254 


Plane 


Y (mm) 


(p (degrees) Q (W) 


A 


7.327 


46 


1.20 


B 


21.370 


58 


6.392 


C 


39 . 688 


70 


15.349 


D 


62.889 


80 


18.388 


E 


87.313 


86 


27 . 628 



108 



OF REFERZNCES 



1. Eilder, D.S., En.-crai”.sr.~ Mcdelinc of Buoyant Mor-enturri 
Jets in Water , Mechanical Engineer's Thesis, Maval 
Postgraduate School, Monterey, CA, 1981. 

2. Nickoden, Deterriination of Velocity Profiles in a 

- Turbulent Jet in Crossflow Using a laser Doppler 

Velocir.eter , M.S. Thesis, Naval Postgraduate School, 
Monterey, CA, 1984. 

3. Oak Ridge Nauional Laboratory, Oak Ridge, Tennes 

Report ORNL-4685 , Analysis of Round, Turbulent, , 

Jets Discharged to Flowing Stratified Arnients , by 
E.A. Hirst. 

4. Incropera, F.P. and Dewitt, D.?., Fundamentals of Heat 
Transfer , John Wiley & Sons, 1981. 

5. Beckwith, T.G., Buck, L.W., and Marangoni, R.D., 
Mechanical Measurements , 2nd ed. , Addison-Wes ley Pub- 
lishing Comtpany , 1971. 

6. Merzkirch, W. , Flow Visualization, Academic Press, 
1974. 

7. Caceci, M.S. and Cacheris, W.p., "Fitting Curves to 
Data," 3vte : The Small Systems Journal, Vol . 9, 

pp. 340-362, May 1984. 

8. Integrated Software System.s Corporation, DISSPLA(R) 
Users Manual , 1981. 

9. Holman, J.P., Exper im.ental Methods for Engineers , 
McGraw-Hill Book Ccmipany , 19 78. 



109 



UJ w 



INITIAL DISTRIBUTION LIST 



No. Copies 
2 



1. Defense Technical Information Center 
Cameron Station 
Alexandria, Virginia 22314 

2. Library, Code 0142 2 

Naval Postgraduate School 

Monterey, California 93943 

3. Department Chairman, Code 69MX 1 

Department of Mechanical Engineering 

Naval Postgraduate School 
Monterey, California 93943 

4. Dr. William G. Culbreth, Code 69Cb 5 

Department of Mechanical Engineering 

Naval Postgraduate School 
Montery , California 93943 

5. Ronald J. Matoushek 5 

R. D. #1 

Waymart, Pennsylvania 18472 



110 



Thesis 

M38075 

c.l 



Thesis 

M38075 

c.l 



O -I 



0384 



Matoushek 

Design and construc- 
tion of a computer 
controlled microtherm- 
couple prohe for the 
study of houyant jets. 



A -4 ^ o » 



Matoushek 

Design and construc- 
tion of a computer 
controlled microtherm- 
couple prohe for the 
study of houyant jets. 



thesM38075 

Design and construction of a computer co 




3 2768 002 12518 9 
DUDLEY KNOX LIBRARY 












